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1.  DESCRIPTION  OF  RESEARCH  AND  RESULTS 


The  research  goal  of  this  project  may  be  summarized  into  the 
following  three  parts 

a)  Waves  forced  by  condensation  heating  in  the  troposphere 

We  have  studies  the  problem  of  the  vertical  structure  of  trop 
ical  waves  generated  by  cumulus  heating  for  two  cases:  l)heating 
is  controlled  by  large-scale  low-level  convergence  (CISK  model) 
and  2)  heating  receives  no  feedback  from  large-scale  motion 
Cforced  model) . In  the  first  case  the  internal  wave-CISK  system 
is  treated  and  upper  and  lower  bounds  for  the  vertical  wavelength 
of  the  unstable  waves  under  normal  heating  conditions  are  esta- 
blished through  analysis  of  the  frequency  equation.  The  lower 
bound  excludes  the  possibility  of  excitation  or  maintenance  of 
short  vertical  wavelengths  relative  to  the  vertical  scale  of  heat 
ing.  The  vertical  structure  of  the  most  unstable  waves  is  also 
computed.  These  results  were  reported  in  Chang  Cl976a) . In  the 
second  case  the  problem  of  scale-selection  of  Kelvin  waves  in  the 
stratosphere  is  analyzed  using  an  inviscid  linear  model.  The  re- 
sults modify  Holton's  (1973)  theory  in  that  1)  the  forcing  is 
most  efficient  for  the  longest  zonal  wavelength  even  if  the  heat 
sources  are  distributed  randomly,  and  2)  the  most  favored  verti- 
cal wavelength  of  the  excited  waves  is  about  twice  the  vertical 
scale  of  heating.  These  results  were  reported  in  Chang  (1976b). 
We  have  also  treated  the  linear  model  with  friction  added  which 
is  called  for  by  several  budget  studies  of  cumulus -related  tropi- 
cal weather  systems,  including  the  monsoon  (Holton  and  Colton, 
1972),  the  cloud  clusters  (Williams  and  Gray,  1973),  and  the 
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easterly  waves  (Reed  and  Recker,  1971;  Wallace,  1971,  Reed  and 
Johnson,  1974;  Chu,  1976).  This  friction  is  believed  to  be  due 
to  vertical  cumulus  transports.  Our  analysis  shows  that  this 
damping  has  a strong  influence  at  low  frequencies  on  the  forced 
equatorial  waves  and  results  in  two  types  of  dispersive  relation- 
ships. The  first  type  is  characteristic  of  the  regular  internal 
gravity  waves  which  have  fast  phase  speeds  and  weak  vertical  atten- 
uation. The  second  type  is  dominated  by  the  viscous  damping  time 
scale  and  has  slow  phase  speed  and  strong  vertical  trapping.  This 
theory  predicts  that  the  stratospheric  oscillations  may  be  identi- 
fied with  the  first  type  and  the  tropospheric  oscillations  with 
the  second  type.  In  the  case  of  Kelvin  waves  the  results  can  be 
used  to  explain  consistently  both  the  observed  stratospheric  Kelvin 
waves  and  the  planetary-scale  Kelvin-like  oscillations  in  the 
troposphere  including  the  40-50  day  oscillation  and  the  winter  mon- 
soon and  Walker  circulations.  These  results  are  described  in 
Chang  Cl977a). 

b)  The  local  barotropic  instability  produced  by  the  spatially 
varying  basic  flow. 

Here  we  have  studied  the  behavior  of  waves  superimposed  upon 
a barotropically  unstable  mean  wind  field,  which  varies  in  x and 
y.  This  mean  wind  field  roughly  simulates  the  200  mb  easterly  jet 
south  of  the  Tibetan  high  during  the  northern  hemisphere  summer. 

The  barotropic  vorticity  equation  is  linearized  and  a linear  damp- 
ing term  is  added.  The  following  mean  u component  is  used 

uCx,y)  = -UCx)  sech^Cy/dCx) ) - , 
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where  dCx)  is  a measure  of  the  width  of  the  jet.  The  quantity 
U(x)  and  u are  chosen  in  such  a way  that  the  mean  flow  is  non- 
divergent.  A rectangular  domain  is  employed  with  a time  periodic 
boundary  contition  on  the  east  and  a radiation  condition  on  the 
west,  which  simulates  the  propagation  of  small  amplitude  waves 
through  the  easterly  jet  region.  The  vorticity  equation  is  solved 
with  the  use  of  finite  differences,  and  when  the  boundary  condi- 
tions are  properly  adjusted  the  waves  move  smoothly  across  the  re- 
gion and  out  the  westr'-n  boundary.  After  a certain  time  interval 
of  numerical  integration,  the  solution  becomes  periodic  everywhere, 
with  the  frequency  that  is  specified  on  the  eastern  boundary.  As 
the  waves  move  through  the  domain  they  grow  or  decay  spatially  in 
reaction  to  the  local  basic  flow,  even  though  at  each  point  the 
variation  is  purely  harmonic.  We  have  developed  a simple  mechan- 
istic analytical  model  which  describes  the  main  features  of  the 
spatial  variation.  The  simple  model  uses  growth  rates  and  phase 
velocities  which  are  computed  from  a parallel  flow  model  with  the 
wavelengths  taken  from  the  full  numerical  model.  These  wave- 
lengths increase  to  a raixiraum  at  the  point  where  the  jet  is  most 
unstable,  and  then  decrease;  the  phase  velocity  has  a similar  be- 
havior. The  maximum  amplitude  in  the  numerical  solution  occurs 
downstream  from  the  most  unstable  part  of  the  jet,  at  the  point 
where  the  growth  rate  in  the  parallel  flow  model  becomes  zero- 
The  above  results  are  grossly  similar  to  that  expected  from  the 
parallel  flow  theory  of  barotropic  instability.  However,  in  the 
unstable  region  the  resultant  structure  of  the  waves  causes  a spa- 
tial growth  rate  greater  than  predicted  by  the  local  growth  rates 
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computed  with  a parallel  flow  model.  In  the  stable  region,  the 
structure  leads  to  a strong  dynamic  damping  which  is  apparently 
due  to  the  presence  of  the  continuous  spectrum  modes  (Case,  1960; 
Pedlosky,  1964).  When  a uniform  advective  velocity  is  added  to  a 
variable  mean  flow,  the  difference  between  the  behavior  of  the 
computed  waves  and  that  implied  by  the  parallel  flow  theory  is 
somewhat  reduced.  However,  in  this  case  a stronger  zonal  asymmetry 
in  the  spatial  growth  rate  with  respect  to  the  jet  maximum  occurs 
as  a result  of  slower  adjustment  of  the  wave  structure  to  the  local 
stability  conditions.  Chang  and  Williams  C1976)  presented  the 
analytic  model  and  discussed  some  other  analytical  aspects  of  the 
problem.  The  numerical  model  results  and  comparisons  with  the 
analytic  model  will  be  reported  by  Tupaz,  Williams  and  Chang  C197S) . 

c)  The  numerical  modeling  of  the  wave-basic  flow  interacting 
system. 

In  this  part  we  adopted  a simple  numerical  model  formulated 
by  Monaco  and  Williams  C1975)  to  simulate  the  planetary-scale  mon- 
soon circulation  and  the  interaction  with  synoptic  waves.  In 
order  to  concentrate  on  the  vorticity  dynamics,  the  heating  is 
specified  with  a horizontal  distribution  given  by  the  time-mean 
200-mb  divergence  of  the  1967  summer  analyzed  by  Krishnamurti  and 
Rogers  (1970).  The  domain  of  the  model  extends  from  18°S  to 
46°N  and  there  are  three  a-levels  in  the  vertical.  Topography 
is  not  included  although  its  thermal  effect  is  implied  in  the 
heating  function.  Two  preliminary  experiments  have  been  carried 
out  to  study  the  vorticity  budget  of  the  planetary  scales  and  the 
non-linear  barotropic  energy  conversions.  In  one  experiment  the 
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heating  is  fixed  in  time  while  in  the  other  it  fluctuates  with  a 
period  of  10  days.  The  time-mean  fields  with  moderate  damping  in 
both  experiments  show  a significant  westward  phase  shift  of  the 
Tibetan  high,  which  suggests  that  the  correlation  of  the  fluctua- 
tions of  the  divergence  and  vorticity  as  suggested  by  Holton  and 
Colton  (1972)  is  probably  not  sufficient  for  the  required  planetary- 
scale  vorticity  sink.  However,  the  integration  period  is  only 
three  times  the  fluctuation  period  in  the  second  experiment  and 
the  amplitude  of  fluctuation  is  about  half  that  observed.  So  the 
full  effect  of  the  fluctuation  may  not  have  been  properly  repre- 
sented in  the  experiment.  A longer  integration  with  larger  fluc- 
tuation is  probably  necessary  to  assess  fully  the  effect  of  the 
fluctuating  heating.  On  the  other  hand,  the  positions  of  the 
simulated  tropical  upper  tropospheric  troughs  in  our  model  agree 
quite  well  with  those  observed.  This  may  suggest  that  damping  due 
to  cumulus  transport  is  more  likely  to  account  for  the  vorticity 
sink  in  the  Tibetan  high  area,  ’Ve  also  found  singificant  transient 
wave  activities  at  the  upper  levels  in  both  experiments  and  they 
appear  to  he  drawing  energy  from  the  planetary  scale  barotropi- 
cally,  in  accordance  with  the  observation  by  Krishnamurti  C.1971)  . 

It  is  possible  that  they  play  an  important  role  in  reducing  the 
amplitudes  of  the  long  waves  which  agree  quite  well  with  those 
observed  despite  the  phase  problem.  This  point  is  borne  out  by 
comparison  with  the  excessively  large  amplitude  found  in  Holton 
and  Colton's  linear  model  when  strong  damping  is  excluded.  In 
addition,  these  short  waves  are  most  active  in  the  fluctuating 
heating  case,  suggesting  that  local  barotropic  instability  at 
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upper  levels  may  be  more  important  than  that  implied  by  the  time- 
mean  wind.  These  results  have  been  reported  by  Chang  and  Williams 
(1976b)  and  Chang  and  Pentimonti  (1977) , 

Some  of  the  results  concerning  the  planetary-scale  monsoon, 
including  the  interpretation  of  the  stationary  equatorial  winter 
monsoon  circulation  as  a viscous,  thermally- forced  Kelvin  wave 
(part  a) , the  preliminary  result  of  the  study  of  barotropic  in- 
stability of  the  upper- tropospheric  summer  monsoon  basic  flow 
(part  b) , and  the  preliminary  result  of  the  planetary  scale  vor- 
ticity  budget  using  forced  numerical  model  Cpai't  c)  are  included 
in  a review  paper  by  Chang  (1977b).  The  research  on  these  prob- 
lems is  being  continued  under  NSF  Grant  ATM77-14821  which  began 
on  15  September  1977. 


REFERENCES 


Case,  K.M. , 1960:  Stability  of  inviscid  plane  Couette  flow.  Phys 
of  Fluids,  2»  143-157. 

Chang,  C.-P.,  1976a:  Vertical  structure  of  tropical  waves  main- 
tained by  internally- induced  cumulus  heating.  J.  Atmos.  Sci., 
32,  729-739. 

Chang,  C.-P.,  1976b:  Forcing  of  stratopsheric  Kelvin  waves  by 
tropospheric  heat  sources.  J.  Atmos . Sci . , 33 , 740-744. 

Chang,  C.-P.,  1977a:  Viscous  internal  gravity  waves  and  low- 

frequency  oscillations  in  the  tropics.  Submitted  to  J.  Atmos. 
Sci . 

Chang,  C.-P.,  1977b:  Some  theoretical  problems  of  the  planetary- 
scale  monsoons.  Pure  and  Appl.  Geophys . , 115  , 1089-1109. 

Chang,  C.-P.  and  R.T.  Williams,  1976a:  Barotropic  instability  of 
a spatially- varying  mean  flow.  Conference  on  atmospheric 
waves  and  stability,  Seattle,  March-April  1976.  Bull.  Amer. 
Meteor.  Soc.,  57 , 101. 

Chang,  C.-P.  and  R.T.  Williams,  1976b:  Study  of  planetary  and 
synoptic-scale  waves  of  the  northern  summer  monsoon.  CARP 
Study  Conference  on  tropical  numerical  modeling,  Exeter, 
England,  April  1976. 

Chu,  J.-H.,  1976:  Vorticity  in  maritime  cumulus  clouds  and  its 

effects  on  the  large-scale  budget  of  vorticity  in  the  tropics. 
Ph.D.  thesis,  UCLA,  123  pp. 

Holton,  J.R.,  1973:  On  the  frequency  distribution  of  atmospheric 
Kelvin  waves.  J . Atmos . Sci . , 30  , 499-501. 

Holton,  J.R.  and  D.E.  Colton,  1972:  A diagnostic  study  of  the 

vorticity  balance  of  200  mb  in  the  tropics  during  the  northern 
summer.  J.  Atmos.  Sci. , 2^,  1124-1128. 

Krishnamurti , T.N.,  1971:  Tropical  east-west  circulations  during 
the  northern  summer.  J.  Atmos . Sci . , 28  , 1342-1347. 

Krishnamurti,  T.N.  and  3.  Rogers,  1970:  200  millibar  wind  field 

June,  July  and  August,  1967.  Rept  70-2,  Dept,  of  Meteorology, 
Florida  State  University,  Tallahassee. 

Monaco,  A.V.  and  R.T.  Williams,  1975:  An  atmospheric  global  pre- 
diction model  using  a modified  Arakawa  differencing  scheme. 
Naval  Postgraduate  School  Technical  Rept.  NPS- 51Nu75041 , 86  pp 

Pedlosky,  J. , 1964:  An  initial  value  problem  in  the  theory  of 
barotropic  instability.  Tellus , 16  , 12-17. 


Reed,  R.  J.  and  R.  H.  Johnson,  1974:  The  vorticity  budget  of  sy- 
notpic-scale  wave  disturbances  in  the  tropical  western  Pacific, 
J.  Atmos.  Sci.,  31,  1784-1790. 

Reed,  R.J.  and  E.E.  Recker,  1971:  Structure  and  properties  of  sy- 
noptic-scale wave  disturbances  in  the  equatorial  western  Pacific 
J.  Atmos . Sci.  , 28 , 1117-1133. 

Tupaz,  J.B.,  R.T.  V/illiams  and  C.-P.  Chang,  1978  : A numerical 

study  of  barotropic  instability  in  a zonally  varying  easterly 
jet.  J.  Atmos . Sci.  ,35 , Cin  press). 

Wallace,  J.M. , 1971:  Spectral  studies  of  tropospheric  wave  dis- 
turbances in  the  tropical  western  Pacific.  Rev.  Geophys . 

Space  Phys . , 9^,  557-612  . 

Williams,  K.T.  and  W.M.  Gray,  1973:  Statistical  analysis  of 

satellite-observed  trade  wind  cloud  clusters  in  the  western 
North  Pacific.  Tellus , 25,  313-336. 


LIST  OF  PUBLICATIONS 


Chang,  C.-P.,  1976a:  Vertical  structure  of  tropical  waves  main- 
tained by  internally- induced  cumulus  heating.  J.  Atmos.  Sci., 
729-739. 

Chang,  C.-P.,  1976b:  Forcing  of  stratospheric  Kelvin  waves  by 
tropospheric  heat  sources.  J.  Atmos . Sci.  , 33 , 740-744. 

Chang,  C.-P.,  1976c:  Comments  on  "Instability  theory  of  large- 
scale  disturbances  in  the  tropics."  J.  Atmos.  Sci., 33, 
1667-1668. 

Chang,  C.-P.,  1977a:  Viscous  internal  gravity  waves  and  low- 

frequency  oscillations  in  the  tropics.  Submitted  to  J.  Atmos. 
Sci . 


Chang,  C.-P.,  1977b:  Some  theoretical  problems  of  the  planetary- 
scale  monsoons.  Pure  and  Appl.  Geophys.,  115 , 1089-1109. 

Chang,  C.-P.  and  R.J.  Pentimonti,  1977:  A numerical  study  of 

time-mean  northern  summer  monsoon  with  steady  and  fluctuating 
heating.  Proc.  Internatl.  Symp.  on  Monsoons,  New  Delhi,  India. 

Tupaz,  J.B.,  R.T.  Williams  and  C.-P.  Chang,  1977:  A numerical 

study  of  the  locally  unstable  barotropic  easterly  jet.  Proc. 
Internatl.  Symp.  on  monsoons.  New  Delhi,  India. 


12 


THESES 


Pentimonti,  R.J.:  Study  of  monsoon  circulation  with  steady  and 
fluctuating  heating.  M.S.  thesis.  December  1976. 

Tupaz , J.B.:  A numerical  study  of  barotropic  instability  of  a 
zonally  varying  easterly  jet.  Ph.D.  thesis,  June  1977. 

SCIENTIFIC  COLLABORATORS 

R.J.  Pentimonti,  graduate  student 

J.B.  Tupaz,  graduate  student 

L.C.  Chou,  graduate  research  associate 

COMMENTS 

This  research  is  being  continued  under  NSF  Grant  ATM77-14821. 

I 


13 


Reprinted  from  Journal  of  the  Atmospheric  Sciences,  V'^ol.  33,  No.  5,  May  1976 
American  Meteorological  Society 
Printed  in  U.  S.  A. 


Vertical  Structure  of  Tropical  Waves  Maintained  by 
IntemaUy-Induced  Cumulus  Heating 

Chih-Pei  Chang 

Department  of  Meteorology,  Naval  Postgraduate  School,  Monterey,  Calif.  939-10 
(Manuscript  received  25  November  1975,  in  revised  form  22  January  1976) 

ABSTRACT 

Solutions  to  the  wave-CISK  (conditional  instability  of  the  second  kind  with  cumulus  heating  being  in- 
duced by  low-level  internal  wave  convergence)  system  are  obtained  to  study  the  vertical  structure  of  margin- 
ally unstable  waves.  A diabatic  heating  profile  is  specified  that  resembles  those  observed  and  those  theoret- 
ically derived  from  simple  parameterization  schemes.  Upper  and  lower  bounds  for  the  vertical  wavelength 
of  the  unstable  waves  under  normal  heating  conditions  are  established  through  analysis  of  the  frequency 
(subility)  equation.  The  lower  bound  excludes  the  possibility  of  excitation  or  maintenance  of  short  vertical 
wavelengths  (relative  to  the  vertical  scale  of  heating)  by  wave-ClSR.  The  calculated  growth  rates  indicate 
that  this  result  is  basically  insensitive  to  the  vertical  heating  profile.  The  vertical  structure  of  the  most 
unstable  waves  is  also  computed  and  the  possible  roles  played  by  CISK  in  large-scale  tropical  waves  are 
discussed  in  light  of  these  results. 


1.  Introduction 

The  possibility  that  large-scale  waves  with  short 
vertical  wavelengths  may  be  excited  in  the  tropical 
atmosphere  has  been  raised  by  Lindzen  (1967),  who 
showed  that  the  equatorially  trapped,  internal  gravity 
and  internal  Rossby  waves  usually  have  a small 
positive  equivalent  depth  and  therefore  a short  ver- 
tical wavelength.  Holton  (1969,  1972a)  later  pointed 
out  that  such  short  vertical  wavelengths,  if  they  exist 
in  the  troposphere,  would  cause  great  difficulties  for 
numerical  weather  prediction  in  the  equatorial  region 
because  of  the  stringent  requirement  for  high  vertical 
resolution.  Whether  short  vertical-scale  waves  actu- 
ally have  a significant  presence  in  the  tropical  tropo- 
sphere or  not  is  unclear  due  to  inadequate  observations. 
Studies  based  on  conventional  radiosonde  data  (e.g., 
Wallace,  1971;  and  others)  tend  to  suggest  that  most 
of  the  variance  of  the  waves  in  the  tropical  tropo- 
sphere is  associated  with  vertical  scales  comparable 
to  the  scale  height.  However,  Madden  and  Zipser’s 
(1970)  analysis  of  high-resolution  rawinsonde  data 
obtained  during  the  Line  Island  Experiment  showed 
that  a quite  short  vertical  wavelength  (~3  km) 
existed  over  the  central  equatorial  Pacific.  Further 
observational  studies  thus  seem  necessary  to  resolve 
this  problem. 

Condensation  heating  due  to  deep  cumulus  con- 
vection has  been  recognized  as  a primary  energy 
source  for  large-scale  tropical  motions.  Two  possible 
mechanisms  which  may  be  responsible  for  the  excita- 
tion of  Wge-scale  tropical  waves  by  heating  have 


been  investigated.  Tnc  first  is  a one-way  forcing  by 
the  cumulus  scale  which  receives  no  feedback  from 
the  large  scale.  Numerical  study  by  Holton  (1971, 
1972a,  1973)  and  analytical  work  by  Chang  (1976) 
have  shown  that  this  mechanism  excites  tropical 
waves  that  have  vertical  scales  comparable  to  those 
of  the  forcing,  which  are  usually  10-14  km.  The 
problem  of  short  vertical  wavelengths  thus  appears 
unimportant  for  forced  waves.  The  second  possible 
mechanism  is  the  conditional  instability  of  the  second 
kind  (CISK)  which  involves  a feedback  to  cumulus 
convection  by  low-level  convergence  of  large-scale 
motion.  This  low-level  convergence  may  be  due  to 
frictionally  induced  cross-isobaric  flow,  as  in  the  case 
of  typhoons  (which  shall  be  called  “Ekman-CISK”), 
or  due  to  the  structure  of  the  internal  waves  (which 
shall  be  called  “wave-CISK”).  The  Ekman-CISK  has 
been  shown  by  Chang  (1971)  and  Chang  and  Piwowar 
(1974)  to  possess  no  preferred  horizontal  scale  for 
synoptic  wave  disturbances.  Thus  Ekman-CISK  may 
not  explain  the  generation  of  synoptic-scale  waves, 
but  it  may  be  quite  important  in  maintaining  the 
waves  once  they  are  initiated.  The  wave-CISK  mecha- 
nism has  been  investigated  by  Hayashi  (1970),  Lindzen 
(1974)  and  Lindzen  el  al.  (1975).  Although  Lindzen 
(1974)  and  Lindzen  el  al.  (1975)  found  no  direct 
horizontal  scale  selection,  the  most  unstable  wave- 
CISK  mode  has  a vertical  scale  seemingly  comparable 
to  the  depth  of  the  sub-cloud  layer  or  the  low-level 
convergent  moist  layer,  depending  on  the  particular 
model  aspects.  A short  vertical  scale  may  therefore 
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be  important  for  the  excited  waves.  Using  the  cal- 
culated vertical  scale  they  also  speculated  that,  due 
to  the  Doppler-shifting  effect  of  the  variable  basic 
flow  and  the  small  growth  rates,  a period  of  ~5  days 
associated  with  a wavenumber  zero,  mi.xed  Rossby- 
gravity  mode  is  likely  to  dominate  the  tropical  spectra. 
This  periodicity  thus  provides  a basic  frequency  for 
tropical  systems.  Depending  on  the  local  basic  flow, 
it  ma\'  force  resonantly  higher  wavenumbers  that  may 
be  relevant  to  the  observed  waves. 

Since  the  growth  rates  associated  with  wave-CISK 
are  usually  very  small,  it  may  be  expected  that  the 
CISK  waves  have  properties  similar  to  neutrally 
forced  waves.  The  difference  in  the  selection  of  vertical 
scales  between  the  forced  models  and  the  wave-CISK 
models  thus  appears  to  be  worth  further  investigation. 
In  addition  to  its  importance  for  tropical  numerical 
weather  prediction,  the  vertical  scale  also  determines 
other  properties  of  the  e.xcited  internal  waves.  The 
purpose  of  this  paper  is  to  examine  the  growth  rates 
associated  with  wave-CISK,  particularlv  with  regard 
to  the  vertical  scale  selection  mechanism. 


vertical  structure  equations  by  assuming  that 
iv  = Y.  Y „iy)wn(,z)  exp[z/ i2H)2, 

n 


where  the  density  factor  exp  [z/(2//)]  is  separated 
for  convenience.  The  meridional  structure  equation  is 


(Py  „ /kii  up  S‘y‘\ 

+( k-^+ ) 

df  \(J1  ghn  ghj 


F„=0, 


(6) 


where  g is  the  gravitational  constant  and  h^,  the 
equivalent  depth,  is  the  separation  constant.  Matsuno 
(1966)  and  Lindzen  (1967)  have  shown  that  the  solu- 
tions to  (6)  which  satisf>  the  boundary  condition 

Fn^  as  lyj— ►»,  (7) 

and  lead  to  the  frequency  equation 
fkH  U^\(^*n)' 

( = 2«+l,  n=-l,0,l,2 (8) 

\ to  ghJ  0 


2.  Basic  equations 

Neglecting  the  shear  of  the  mean  zonal  wind,  the 
linearized  zonal  and  meridional  momentum,  hydro- 
static, thermodynamic  energy  and  continuity  equa- 
tions for  a single  zonal  wavenumber  k on  an  equatorial 
beta-plane  may  be  written 


Y.{y)  = 


-//„_,(i)— 


1— (IA,)‘ 

to 


r 1 

2 l+-{gh„)i 

- to 


-^,+l({) 


Xexp(-{V2), 


iu:il—0yv=  ~ik0, 

(1) 

(2) 

iuiip-0\u  , 

dy 

Rt 

— 1 

(3) 

dz  H 

ia)f-f-wr=— , 

(4) 

di  d 

ika^ — («-*"'u!)=o. 

(S) 

dy  dz 

where  H,  H,  tb,  T,  0 and  (5  are  the  perturbation  zonal 
velocity,  meridional  velocity,  vertical  velocity,  tem- 
jjerature,  geopotential  and  diabatic  heating  rate, 
respectively;  to  is  the  Doppler-shifted  frequency, 
H a constant  scale  height,  F the  static  stability, 
Cf  the  specific  heat  at  constant  pressure,  R the  gas 
constant,  0 the  meridional  gradient  of  the  vertical 
compionent  of  earth’s  vorticitv,  y the  meridional  co- 
ordinate, and  z=—H\n  (p/po)  is  the  vertical  coordi- 
nate with  p the  pressure  and  pn  a reference  pressure. 

Eqs.  (l)-(5)  may  be  combined  into  a single  equa- 
tion in  vb  that  may  be  separated  into  meridional  and 


where  (=0*(g^in)~h  and  the  Hermite  polynomials 
have  values  only  for  «^0.  The  meridional 
velocity  solution  is  i',(y) «/7„({)  exp  ( — J’/2).  Lamb 
(1973)  has  solved  (6)  using  the  more  general  form 
of  the  confluent  hypergeometric  functions  without  the 
restriction  of  (7).  In  such  a case,  the  lateral  boundary 
conditions  are  given  by  the  choice  of  vanishing  me- 
ridional velocity  at  some  finite  |yl  where  the  con- 
fluent hypergeometric  functions  have  a node.  The 
Hermite  polynomials  are  a sub-set  of  the  confluent 
hyjjergeometric  functions. 

The  vertical  structure  equation  which  is  most  rele- 
vant in  the  present  calculation  is 


where 


dhVn  Qn 

|-X»®U/,= , 

dz*  ghn 

S 1 

ghn  4H* 


(9) 

(10) 
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— r, 
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(5,F,-'exp[-z/(2H)]. 
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The  parameter  A,  is  a measure  of  the  vertical  wave- 
number  and  is  generally  complex  for  unstable  waves. 
In  the  case  of  forced  modes  w and  k are  known  pa- 
rameters of  the  forcing  function  so  that  //»  and  A,  are 
obtained  from  (8)  and  (10)  respectively.  For  CISK 
modes,  on  the  other  hand,  A„  is  determined  as  an 
eigenvalue  of  (9)  and  w is  then  found  from  (8)  and  (10). 
In  the  present  calculation,  the  heating  function  Q'  is 
assumed  to  possess  a white-noise  distribution  for  all 
Hermite  modes  (all  values  of  n)  so  that  the  calcula- 
tion of  h will  be  independent  of  n.  The  subscript  n 
will  therefore  be  dropped  in  the  subsequent  discus- 
sions. It  is  well  known  that  the  vertical  structure 
equation  (9)  can  also  be  derived  by  considering  two- 
dimensional  wave  motion  only. 

As  in  the  forced  model  by  Chang  (1976)  the  fol- 
lowing boundary  conditions  are  used  to  solve  (9) ; 


K = 


0, 


at  z = 0 
at  3=z, 


(lla) 

(llb) 


where 


Cl  — fC2. 


Here  z,  is  the  height  of  tropopause,  Ci,  C.  are  con- 
stants, and  the  condition  (11b)  results  from  the  re- 
quirement that  latent  heating  vanishes  at  z,.  The 
parameter  r is  a reflection  coefficient  which,  in  the 
absence  of  vertical  wind  shear,  is  given  by  the  speci- 
fication of  the  static  stability  distribution.  If  the 
heating  is  assumed  to  vanish  below  the  cloud  base  (z^), 
the  solution  to  (9)  may  be  written  in  the  form 


r"  Q' 

a'(z)  = / sinAz — i/z,  z^Zi,  (12a) 

XCl-fr)  J..  gh 


sinAz  /•”  Q' 

a(z)  = / (re***-|-e~‘^*) — dz 

A(l-f-r)y,  gh  A(l-hr) 


r . 

C / smAz — dz, 

J.. 


z,>z>Zc,  (12b) 


sinAz  /•'•  O' 

w(z)  / (re’^'-f  «-•»')— rfz, 

gh 


Zc^z^O.  (12c) 


The  solution  procedure  using  the  Green’s  function 
technique  is  the  same  as  in  Lindzen  (1974). 


3.  The  heating  function 

The  parameterization  of  the  effect  of  condensation 
heating  due  to  cumulus  convection  is  a major  problem 
in  tropical  modeling.  Important  contributions  on  this 
problem  have  been  made  by  the  sophisticated  schemes 
develofjed  by  Arakawa  and  Schubert  (1974),  Ooyama 
(1971),  and  others.  On  the  other  hand,  in  simplified 


analytical  models,  simpler  schemes  must  be  used  to 
keep  the  mathematical  analysis  tractable.  Usually  in 
these  models  the  vertical  distribution  of  the  heating 
is  specified  as  a given  function  that  is  invariant  in 
time.  The  major  criticism  of  this  type  of  parameteriza- 
tion is  that  the  results  are  sensitive  to  the  specified 
heating  profile  such  that  one  can  obtain  the  desired 
result  by  simply  tuning  the  heating  profile.  Another 
criticism  is  that  the  inherent  time-variation  of  the 
heating  function  due  to  the  complicated  scale  inter- 
action cannot  be  properly  includ^.  However,  recent 
observational  studies,  using  both  the  spectral  analysis 
technique  and  the  composite  technique,  have  produced 
a consistent  picture  of  the  vertical  heating  profile  due 
to  latent  heat  release  over  the  tropical  oceanic  areas. 
The  profiles  observed  by  Nitta  (1970),  Wallace  (1971), 
Reed  and  Recker  (1971)  and  Williams  and  Gray  (1973) 
are  plotted  in  Fig.  1.  It  can  be  seen  that  in  each  case 
a maximum  occurs  in  the  middle  troposphere  between 
the  6 and  9 km  levels.  On  the  other  hand,  observa- 
tional analysis  of  the  divergence  and  vertical  velocity 
by  Schubert  and  Reed  (1975)  using  recent  GATE 
data  has  yielded  a result  that  implies  a vertical  heating 
profile  different  from  those  shown  in  Fig.  I.  They 
found  that  the  maximum  vertical  velocity,  which 
usually  coincides  with  the  maximum  heating,  based 
on  scaling  arguments  (Holton,  1972b),  occurs  at  a 
much  lower  level.  However,  analyses  of  the  same  data 
set  by  Petrossiants  et  al.  (1975)  and  .Antsipovich  et  al. 
(1975)  of  the  USSR  Hydrometeorological  Center  have 
resulted  in  profiles  which  are  in  close  agreement  with 
Fig.  1.  In  any  case,  the  profiles  included  in  Fig.  1 are 


Fig.  1.  Vertical  heatinK  functions  deduced  from  observations  and 
from  simple  theoretical  parameterization  schemes. 
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based  on  a much  larger  sample  of  data  and  may  be 
viewed  as  an  indication  that  the  tropical  atmosphere, 
on  a statistical  basis,  tends  to  maintain  a similar 
vertical  profile  for  the  large-scale  heating  function. 

Results  of  a theoretical  calculation  by  Kuo  (1965) 
using  his  parameterization  scheme  are  also  included 
in  Fig.  1.  Kuo’s  parameterization  was  originally  based 
on  the  assumption  of  horizontal  mixing  between  cloud 
air  and  environmental  air.  Recently  Kuo  (1974)  has 
shown  that  his  formulation  is  equivalent  to  the  con- 
sideration that  heating  is  due  to  the  adiabatic  com- 
pression of  the  descending  environment.  Kuo’s  scheme 
has  been  used  in  many  tropical  numerical  models 
with  a degree  of  success  not  inferior  to  any  other 
parameterization  scheme.  It  is  therefore  interesting  to 
compare  his  heating  profiles  with  those  observed.  In 
Fig.  1,  the  Kuo-1  curve  is  for  the  case  of  no  entrain- 
ment and  the  Kuo-2  curve  is  for  the  case  of  moderate 
entrainment.  It  is  clear  that  the  latter  curve,  having 
a maximum  near  7 km,  agrees  quite  well  with  those 
deduced  from  observations.  The  average  healing  profile 
obtained  from  the  GFDL  general  circulation  model 
as  calculated  by  Hayashi  (1973)  is  also  included, 
because  the  GFDL  model  uses  a different  parame- 
terization scheme:  the  moist  convective  ac.'ustment. 
Hayashi  decomposed  the  heating  function  into  com- 
ponents of  various  equatorial  beta-plane  modes.  .Among 
them,  the  profile  for  the  Kelvin  waves  («=  — 1)  at 
the  equator  is  one  of  the  curves  least  resembling 
those  observed.  This  profile  is  also  plotted  in  Fig.  1. 
It  is  seen  that  even  for  this  highly  simplified  parame- 
terization scheme  the  maximum  heating  still  occurs 
near  6 km,  within  the  range  of  the  other  curves. 


HEATING  FtmCTION  0(2) 


21 


Q'di  — 

Fio.  2.  The  specified  heatinir  function  Q’(s). 


Based  on  the  above  discussion  the  heating  function 
for  the  wave-CISK  parameterization  is  specified  as 

()'(z)  = iwA'5,W6e°'' sin  its', 

^'=0,  z>z,  or  z<Zt 

where 


and  Az  = z,— Zf  is  a measure  of  the  vertical  scale  of 
the  heating.  The  coefficient  m specifies  the  strength 
of  the  heating  and  the  factor  J is  needed  because 
the  amplitude  of  the  Fourier  component  is  one-half 
of  the  heating  maximum,  as  only  positive  conden- 
sation heating  is  permitted.  The  heating  is  propor- 
tional to  the  moisture  convergence  in  the  moist  layer 
which  is  represented  by  the  vertical  velocity  Wb  at  the 
top  of  the  moist  layer.  The  tropospheric  value  of  the 
static  stability  St  is  included  in  the  proportionality 
constant.  The  parameter  a is  used  to  vary  the  ma.xi- 
niuni  heating  level  to  test  the  sensitivit.v  of  the  model. 
The  coefficient  .V  is  a normalization  factor  so  that 
when  other  conditions  are  equal  the  total  amount  of 
heat  release  weighted  by  density  in  a column  would 
remain  the  same  with  different  values  of  a.  Thus, 


f e.\pC-3/(2 
J 0 


:/(2fl)],/z' 


/ e.\p[az'— z/  (2//)]  simrz'i/z' 

Jo 

2//{  1 — exp[— Az/  (2tf )])  { [a— Az/ (2//)]=-|-tr=} 
IT  { exp[  j — Az/ (2/f ) ] -t- 1 ) 


Profiles  of  (13)  multiplied  by  the  density  factor 
exp  (z/2//)  are  plotted  in  Fig.  2 for  various  values 
of  a.  The  range  of  — l^u^O,  with  maximum  heating 
between  7 and  9 km,  appears  to  be  representative  of 
the  observed  and  theoretical  profiles  shown  in  Fig.  1. 


4.  Analysis  of  the  stability  equation 

If  the  mixed  layer  with  top  at  Zt  (<z.)  is  assumed 
to  be  the  moisture  layer  which  provides  the  con- 
vergence for  CISK,  the  most  relevant  solution  in 
studying  the  instability  is  (12c)  evaluated  at  Zb.  The 
solution  is 


sinXzt  mwb  St  tAz 

Wb= 

X l+r  gh  2 

re- (e« -|- 1 ) -f  re**  1 ) 1 

=<[ ] 

where  m-mX  and  q = a—i\Az.  Equating  the  coeffi- 
cients of  Wb  on  both  sides  and  using  (10)  leads  to 
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the  stability  equation 


(x^+— )l 
V 4//V 
X 


(!+'•)  (9'^+ir-) 


-+1=0, 

(14) 


gives  two  equations : 

Re(^')=0,  (18) 

Im  ('{'')  =0. 

For  <J  = 0,  (18)  may  be  expanded  to  give 


Re('I'')  = fxr^-Xi''H Vexp(\iZ,)cos,\,s, 

\ 4/7V 

+ exp(\,Sr)cOSX,2r]  + 2XrXiCeXp(X,Z,)sin,\rZ, 


where  we  have  assumed  that  (sin  Xsfc)/X«Z4.  It  can 
be  shown  that  ^^+^“=0  is  not  a solution  to  (l4). 
Eq.  (14)  must  be  solved  for  the  complex  eigenvalue  X 
using  either  numerical  or  graphical  methods.  However, 
before  doing  so  some  valuable  information  regarding 
the  range  of  vertical  scale  of  the  unstable  solutions 
may  be  obtained  by  analyzing  (14). 

The  complex  eigenvalue  X may  be  written  as  X=Xr 
+!X(,  where  the  real  part  X,  is  the  vertical  wave- 
number  while  the  imaginary  part  X,  gives  the  growth 
rate  of  the  waves.  If  the  positive  X,  is  chosen,  the  Ca 
component  in  (11b)  represents  the  downward  phase 
or  upward  energv  propagation  and  the  Ci  component 
represents  the  upward  phase  or  downward  energy 
propagation.  It  then  follows  from  (10)  and  (8)  that 
if  Xi  is  negative,  the  imaginary  part  of  the  frequency  u 
will  be  negative  and  the  waves  will  grow  e.xponentially 
in  time.  Based  on  observations,  and  justified  a pos- 
teriori, we  shall  consider  only  the  case  of  marginally 
unstable  waves,  or  very  small  growth  rates.  Thus  the 
conditions  applied  to  (14)  may  be  written 

Xi=-lX.-l,  (15) 

X.*«X,^  (16) 

In  addition,  we  shall  assume 

X,*»— , (17) 

4/72 

which  is  valid  for  all  practical  purposes  since  the 
scale  height  H of  the  atmosphere  is  ~7  km  and  (17) 
only  places  a restriction  of  ~88  km  for  the  maximum 
vertical  wavelength. 

Eq.  (14)  will  now  be  analyzed  by  assuming  that 
<i=0  and  that  the  stratosphere  {z>z,)  and  troposphere 
(z^z,)  have  the  same  constant  static  stability  S,. 

In  this  case  r = 0 and  since  the  stability 

equation  may  be  re-written  as 

'k'  f X)  = 'k  (X)  (v^+T^)  = -mzb^z(\^+—\ 

2 \ 4//V 

X exp  ( — /X,2«) -t- (/*-!- = 0. 
.Separating  'k'(X)  into  its  real  and  imaginary  parts 


+e.xp  (X,'2c)sinX,zJ-|- 


2[7r=-A2nx,*-X,=)] 

irmziAz 


= 0. 


Simplifying  this  equation  with  the  conditions  (16) 
and  (17)  and  using  (15)  lead  to 

Re('k') 

» {\?  cosXrZ,  — 2Xr I X, , sinXrS,)e.xp(—  |X,|2c) 

I 


-)-(Xr^  co,sX,2c  — 2X,.!X,|sinXr2,)exp(—  |X,  |2c) 

U 


2(-s^-Az^K^) 

-I- =0.  (19) 

inhZh^z 

III 


Since  JcosX^z,  [,  isinXrS,',  !cosX,z„  and  !sinX,2c! 
are  all  ^1  and  exp(- [X,! 2„)<  1,  from  (16)  we  have, 
allowing  for  small  error  under  extreme  cases,’ 


|I|^X,^ 
illl^X,*  ■ 


(20) 


Thus  for  (19)  to  be  satisfied  we  may  write 

Ma.x[Re(^'')]  = Max(I)+Max(II)+III>0, 

where  Max  denotes  the  maximum  value  using  the 
range  (20) ; 

Max(I)  = Max(II)=Xr». 

Thus  we  have 

2ir  2A2Xr’’ 

2Xr^+ ^>0. 

inzkAz  rmzb 

which  yields 


SziAz—trinZb) 

(21) 

A2>irm24. 

(22) 

'The  precise  maximum  for  terms  I and  II  is  1 1 (l-l-4«’)cos 
[tan-'(-2«)|]  e.xp(-|X.U,)|X,®,  where  X,  X,!  Inequality 
(20)  is  true  when  For  «-0.1  which  is  our  assumption  the 
maximum  is  [1.02  expl  — IX, U,)lX,’.  Obviously  if  'X,  liecomes 
laroe  the  exponential  factor  would  become  much  less  than  1. 


If 
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Fig.  3.  Values  of  Re('t)  and  Inif*)  as  a function  of  X,[k'm  ']  and  The  crosses  in  the  Re(4') 

diagram  are  solutions  for  which  Re('k)  = Im('I')=0. 


Choosing  the  minimum  values  (Min)  in  (20),  we  have 
MinCRe(4'')]  = Min(I)+Min(II)+III<0, 


or 


giving 


2ir  2A3X,* 

-2\r-+ 3— <0, 

ffiziAn  TinSb 


. (23) 

Az{Az+TrmZb) 


so  the  vertical  wavelength,  L = 2ir/\r,  is  limited  by 
the  following  bounds: 

13  km<2)<40  km. 

Thus  waves  with  a vertical  wavelength  shorter  than 
13  km  will  not  be  excited  by  wave-CISK.  If  such 
waves  are  first  excited  by  another  mechanism,  the 
above  result  indicates  that  they  cannot  be  maintained 
by  wave-CISK. 


So  (21)  and  (23)  give  the  upper  and  lower  bounds 
for  \r.  It  is  interesting  to  note  that  X,— ^rr/Az  in  the 
limit  of  m— *0,  corresponding  to  a vertical  wavelength 
of  twice  the  heating  scale,  which  is  the  most  efficiently 
excited  wavelength  for  forced  waves  (Chang,  1976). 
However,  X,  = ir/A3  cannot  be  a solution  for  any 
finite  m because  it  does  not  satisfy  (14). 

Now  we  shall  use  the  values  z,  = IS  km,  Zc=l  km, 
Az=14  km,  z»  = 0.4  km,  w=5.8  which  leads  to  m='9. 
The  value  of  m is  calculated  by  fitting  (13)  to  Reed 
and  Recker’s  (1971)  observed  heating  and  vertical 
velocity  profiles  for  the  wave  trough,  for  which  an 
estimate  of 

»«z»»2.32km,  (24) 

can  be  made  assuming  that  5|=1.22X10~^  s“*  based 
on  a tropospheric  value  of  T — d'C  km“‘.  Eq.  (24)  is 
a good  approximation  as  long  as  Z4<4  km,  so  the 
particular  choice  of  Zb  is  not  crucial.  From  Reed  and 
Recker’s  data  Ws«0.29  cm  s“‘  so  that  the  maximum 
heating  is  ~6.8®C  day“‘.  Substituting  these  values 
into  (21)-(23)  we  obtain 

»•  T 

=0.512>  lx,!  > -0.167  km-‘ 

6.14  18.8 


5.  Eigenvalues  and  eigenfunctions 

a.  Growth  rates 

The  stability  equation  (14)  will  now  be  solved  using 
the  graphical  method  with  the  accuracy  of  the  solu- 
tions improved  by  a numerical  iteration  technique. 
The  graphical  method  is  illustrated  in  Fig.  3,  where 
the  real  and  imaginary  parts  of  'F(X)  given  by  (14) 
with  r=0,  a = 0 and  w = 6 are  plotted  as  a function 
of  X,  and  X,-.  Here  the  constants  used  are  ff  = 7 km, 
Zi=15  km,  Zc  = l km  and  Zb  = 0A  km.  The  intersections 
of  the  zero  Re  (4')  and  zero  Im(4')  lines  give  the 
eigenvalue  solution  X.  It  is  seen  that  there  exist  many 
solutions  with  the  most  (and  only)  unstable  solution 
having  the  smallest  X,  or  the  longest  wavelength.  The 
solutions  are  labeled  by  an  integer  k with  increasing  k 
indicating  increased  X,  or  decreased  vertical  wavelength. 

The  four  longest  vertical  wavelength  solutions  are 
shown  in  Fig.  4 as  a function  of  heating  strength  m. 
It  is  seen  that  |Xj|«Xr  and  that  the  growth  rates  as 
represented  by  — Xi  increase  with  m while  the  vertical 
wavelengths  are  almost  constant.  This  diagram  also 
indicates  that  the  /&=!  solution  is  most  unstable 
throughout  the  range  m = 4 to  m— 10,  which  cor- 
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Fig.  4.  The  four  longest  vertical  wavelength  solutions  as  a 
function  of  the  heating  rate  hi.  The  solid  curves  are  X,  and  the 
dash-dotted  curves  X,.  The  corresponding  vertical  wavelength 
Z.i“2ir/X,  is  given  in  the  parentheses  of  the  X,  coordinate.  The 
ratio  of  static  stabilities  between  stratosphere  and  troposphere 
is  I. 

responds  to  maximum  heating  rates  of  4.7  and  11.7°C 
day“‘  if  «i(,  = ().29  cm  s~*.  In  fact,  it  is  the  only  un- 
stable solution  up  to  tn  = 9.  For  this  solution  the  value 
of  X for  m = 5.8  is  0.388— 0.026i  Ckm“‘3,  which  is  well 
within  the  range  predicted  by  (21)  and  (23). 

The  sensitivity  of  the  eigenvalues  with  respect  to 
the  vertical  heating  profile  (13)  is  shown  in  Fig.  5 in 
which  the  ^ = 1 to  k=4  solutions  with  w = 6 are  given 
as  a function  of  a.  It  is  clear  that  the  vertical  wave- 
lengths are  all  insensitive  to  a.  For  the  growth  rates 
the  k—l  solution  is  almost  independent  of  a,  but  the 
other  solutions  all  increase  with  a.  There  is  an  indica- 
tion that  as  a increased  to  >5  the  shortest  vertical 
wavelength  becomes  the  most  unstable  solution.  How- 
ever, for  the  realistic  range  of  1,  which 

represents  the  statistical  mean  state  of  the  heating 
function  in  the  tropical  atmosphere,  fe  = l is  the  only 
unstable  solution. 

Solutions  for  a different  static  stability  S,  in  the 
stratosphere  are  obtained  by  including  the  radiation 
condition  for  z>3,  as  the  third  boundary  condition. 
In  this  case  the  reflection  coefficient  is  derived  by 
continuity  requirements  of  w and  0 across  z, : 

l-(5./5,)‘ 

r- exp(— 2/X<2,),  for  5,>5,. 

H-(5./5,)» 

The  result  for  5./5i=»2.5  is  shown  in  Fig.  6.  It  is 


Fig.  5.  .Xs  in  Fig.  4 e.’ccept  that  the  solutions  are  given  as  a function 
of  the  heating  profile  parameter  a. 


seen  that  the  general  conclusions  drawn  from  Fig.  5 
remain  valid,  although  waves  with  a shorter  vertical 
wavelength  become  more  unstable  than  the  k=l  waves 
beginning  at  a = 2 which  corresponds  to  a ma.ximum 


Fig.  6.  As  in  Fig.  4 except  that  the  solutions  are  given  as  a 
function  of  the  heating  profile  parameter  a and  that  the  ratio  of 
static  stabilities  between  stratosphere  and  troposphere  is  2.5. 


736 


JOURNAL  OF  THE  ATMOSPHERIC  SCIENCES 


Volume  33 


heating  level  of  ~12  km  or  200  mb,  much  above  all 
observed  levels.  So  our  result  is  basically  insensitive 
to  the  vertical  heating  distribution. 

b.  Vertical  profiles 

The  vertical  structure  of  w,  corresponding  to  the 
X for  1 and  a=0  in  Fig.  6,  is  shown  in  Fig.  7 for  two 
heating  intensities  m=6  and  10.  It  can  be  seen  that 
the  tropospheric  amplitude  reaches  a ma.ximum  be- 
tween 6 and  9 km  and  has  a minimum  at  the  tropo- 
pause  for  both  cases.  In  the  stratosphere  the  amplitude 
for  w = 6 increases  at  a slow  rate  with  height  due  to 
the  density  factor,  but  the  kinetic  energy  actually 
decreases  with  height.  This  decrease  is  due  to  the 
instability  of  the  waves  which  results  in  Xi<0.  The 
decrease  of  the  wave  kinetic  energy  is  much  faster 
for  the  more  unstable  ;«  = 10  case,  so  that  even  its 
amplitude  decreases  with  height.  The  vertical  phase 
distributions  for  both  cases  are  almost  the  same  and 
indicate  that  there  is  virtually  no  vertical  propagation 
in  the  troposphere  in  contrast  to  a downward  propa- 
gation in  the  stratosphere.  The  vertical  wavelength 
computed  from  27r  X,  is  16.2  km  in  the  troposphere 
and  10.2  km  in  the  stratosphere,  because  the  value 


Fig  7,  .Xmplitude  and  phase  of  li  as  a function  of  heii?ht  for 
m — 6 'solid  lines)  and  ih«10  Plotted  lines).  The  phase  lines  for 
lioth  cases  are  approximately  the  same. 


Fig.  8.  .Xs  in  Fig.  7 except  for  t. 


of  X,  is  increased  in  the  stratosphere  by  appro.ximately 
a factor  of  as  can  be  seen  from  (10). 

Vertical  structures  of  other  variables  can  be  com- 
puted utilizing  the  basic  equations  (l)-(5)  and  the 
solutions  to  the  meridional  structure  equation  (6). 
Here  only  the  temperature  structure  for  the  Kelvin 
wave  case  (n=— 1)  will  be  shown.  From  (4)  and  (8) 
we  have 


It  is  obvious  from  (25)  that  the  amplitude  of  the 
temperature  perturbation  decreases  as  the  zonal  wave- 

number  k increases.  The  t„ i structure  corresponding 

to  the  w in  Fig.  7,  is  shown  in  Fig.  8.  It  can  be  seen 
that  the  amplitude  distribution  is  quite  sensitive  to 
the  strength  of  the  heating.  The  discontinuity  at  z,  is 
due  to  the  layered  static  stability  distribution  assumed 
in  our  model.  The  large  temperature  perturbations  in 
the  stratosphere  are,  of  course,  due  to  the  large  static 
stability  there.  Figs.  7 and  8 are  constructed  using 
nondimensional  units.  The  two  figures  may  be  com- 
pared by  letting  the  unit  for  the  w amplitude  be 
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1 cm  s~';  in  this  cast  the  corresponding  unit  for  the 
t amplitude  will  he  [°C],  where  a = kX  (radius  of 
the  earthl  is  the  numher  of  waves  around  the  equator. 
For  these  units  the  amplitude  of  w is  quite  close  to 
the  typically  observed  values  of  "matured”  tropical 
waves,  with  maximum  cm  s“'  for  m = 6 and 
^4.2  cm  s~'  for  w=l().  On  the  other  hand,  the  maxi- 
mum t perturbation  in  the  troposphere  is  almost 
1()°C  for  wavenumber  (a)  1 and  even  larger  in  the 
stratosphere.  Lindzen  (personal  communication)  has 
suggested  that  these  enormously  large  temperature 
fluctuations  indicate  that  large  damping  effect  of 
vertical  momentum  transport  due  to  cumulus  con- 
vection must  be  important.  However,  for  shorter  waves 
the  temperature  perturbation  is  proportionally  smaller. 
If  we  consider  a zonal  wavelength  of  4(X)0  km,  then 
a=l()  and  the  tem])erature  fluctuations  are  all  <1°C 
in  the  troposphere,  in  good  agreement  with  observa- 
tions and  with  Holton  O'^ni’s  numerical  calculations. 
In  fact,  this  magnitude  estimate  may  explain  why 
atmospheric  Kelvin  waves,  which  have  a T amplitude 
~2°C  in  the  lower  stratosphere,  are  not  observed  in 
the  troposphere.  For  such  waves  the  vertical  velocity 
in  the  troposphere  would  be  an  order  of  magnitude 
smaller  than  for  the  easterU'  waves  which  have  a 
typical  zonal  wavelength  km  and  a maximum 

vertical  velocity  ~3  cm  s~'.  On  the  other  hand,  damp- 
ing effect  of  a certain  magnitude  must  play  a role  for 
the  easterl>'  waves  to  explain  the  complete  absence 
of  them  in  the  stratosphere.  Although  these  estimates 
are  based  on  Kelvin  waves  only,  for  other  wave  modes 
a relationship  between  the  tropospheric  temperature 
amplitude  and  the  zonal  scale  must  also  exist  at 
latitudes  somewhat  away  from  the  equator  because 
of  the  following  scaling  arguments.  Wallace  (1971) 
and  Holton  (1972b)  have  shown  that,  for  weak, 
tropical,  synoptic-scale  motions  with  a depth  scale 
comparable  to  the  scale  height,  a diabatic  heating 
rate  of  5-10°C  day~‘  is  almost  completely  balanced 
by  the  adiabatic  cooling.  Thus  the  temperature  fluc- 
tuation is  one  order  of  magnitude  smaller  than  either 
term.  This  balance  may  be  attributed  to  the  smallness 
of  the  Rossby  number  which  is  ~0(1)  or  larger  for 
tropical  synoptic-scale  motions.  For  planetary-scale 
motions  the  length  scale  is  one  order  larger  so  that 
the  Rossby  number  is  ~O(10~’)  at  intermediate 
tropical  latitudes,  and  the  scale  analysis  indicates  that 
the  pressure  and  temperature  perturbations  must  be 
increased  by  one  order  of  magnitude. 

The  phase  diagrams  in  Figs.  7 and  8 indicate  that 
in  the  lower  troposphere  and  the  entire  stratosphere 
the  temperature  leads  the  vertical  velocity  by  ~ i cycle. 
In  the  middle  and  upper  troposphere  (7-12  km  for 
m = 6 and  5-12  km  for  w=10),  they  are  in  phase 
e.xcept  near  the  tropopause  (12-14  km),  where  they 
are  out  of  phase.  These  phase  relationships  are  con- 
sistent with  the  energetics  because  the  conditions  for 


The  latter  condition  implies  that  the  diabatic  heating 
rate  e.xceeds  the  adiabatic  cooling  rate  which,  of 
course,  is  responsible  for  the  continuous  growth  of  the 
waves  with  time. 


c.  Equatorial  trapping  scales 

The  meridional  scale  of  the  unstable  waves  may  be 
estimated  from  the  Hermite  solutions  of  the  meridional 
equation  (6).  The  north-south,  e-folding  width  is 

y,  = Re[2(g/i)Vd]*, 

and  the  ma.ximum  amplitude  for  the  mi.xed  Rossby- 
gra-ity  mode  («  = 0)  occurs  at 

yo=Re[(,?/t)V/3]*. 

The  maximum  for  the  Kelvin  mode  occurs  at  the 
equator.  From  results  obtained  for  the  ^=1  waves, 
/\a:0.388— 0.026/ [km~‘]  so  that  (?/0*“25.5  — 1.7/ [m 
s“‘]  from  (10).  It  follows  that  y,»1480  km  and  yo 
«1050  km  with  /?=■  2.21X  10-“  nf' s-'.  Both  of  the.se 
are  reasonable  values  in  terms  of  the  trapping  of 
waves  on  the  equatorial  beta-plane. 

6.  Concluding  remarks 

We  have  shown  that  tropical  waves  cannot  remain 
unstable  with  respect  to  wave-CISK  if  the  vertical 
wavelength  is  much  smaller  than  the  vertical  scale 
of  heating.  This  result  holds  for  a number  of  reason- 
able heating  profiles  that  have  been  observed  or 
theoretically  derived.  The  dispersive  relationship  of 
the  growth  rates  of  various  equatorial  beta-plane 
modes  has  been  shown  by  Lindzen  (1974)  to  be  un- 
satisfactory in  selecting  the  horizontal  scale  of  tropical 
waves.  From  the  computed  eigenvalues  it  can  be 
shown  that,  even  under  the  extremely  intense  heating 
case  of  m=W,  the  growth  rates  for  Rossby  waves, 
mi.xed  Rossby-gravity  waves,  and  the  longer  zonal- 
scale  Kelvin  waves  are  all  very  small.  This  result  is 
similar  to  that  obtained  by  Chang  and  Piwowar 
(1974)  for  the  Ekman-CISK  and  indicates  that  wave- 
CTSK  is  not  responsible  for  the  initial  selection  mecha- 
nism of  the  tropical  waves.  Once  the  waves  are  excited 
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by  some  other  mechanism,  however,  the  internal  wave 
convergence  could  induce  organized  cumulus  convec- 
tion and  release  wave-CISK  as  an  energy  source.  The 
present  calculation  indicates  that  under  such  cases 
the  shorter  vertical-scale  waves  will  still  not  be  main- 
tained by  wave-CISK.  There  is  no  definite  e.xplana- 
tion  to  the  fact  that  the  gravity  modes  and  the  shorter 
zonal  scale  Kelvin  modes,  which  are  usually  not 
observed  in  the  tropical  atmosphere,  are  most  unstable. 
However,  these  waves  have  ver\-  short  periods  and 
various  arguments  have  been  given  against  their  oc- 
currence on  a regular  basis.  For  example,  if  one  con- 
siders the  lifting  process'*  of  CISK  (Ooyama,  1969; 
Holton,  1972b),  the  short  periods  will  be  insufficient 
for  a low-level  air  parcel  to  reach  its  condensation 
level.  In  any  case,  the  result  that  wave-CISK  can 
neither  excite  nor  maintain  waves  with  short  vertical 
wavelengths  remains  valid  independent  of  these  con- 
siderations. The  concern  about  short  vertical-scale 
excitation  due  to  wave-CISK  in  a tropical  numerical 
model  may  thus  be  relieved.  The  resonant  forcing  of 
tropical  waves  by  a fundamental  wave-CISK  fre- 
quency associated  with  the  zonally-symmetric,  mixed 
Rossbv-gravity  mode  as  hypothesized  b\-  Lindzen 
(1974)  also  does  not  appear  to  be  plausible,  because 
the  equivalent  depth  would  be  too  large  to  give  a 
fundamental  periodicity  near  5 days. 

The  calculated  vertical  structure  of  w and  T applies 
to  all  waves  maintained  by  heating  and  is  not  re- 
stricted to  wave-CISK.  The  eigenvalue  \ found  by 
the  wave-CISK  equation  used  here  may  alternatively 
be  viewed  as  being  excited  by  a specilied  forcing 
[as  in  the  model  by  Chang  (1976)]  with  a small 
linear  damping.  The  damping  coefficient  would  then 
appear  in  the  complex  Doppler-shifted  frequency  in 
the  same  way  an  an  imaginary  frequency  due  to 
instability.  The  computed  magnitudes  of  w and  T 
suggest  that  the  observed  planetary-scale  Kelvin 
waves  are  likely  to  be  a result  of  forcing  rather  than 
instability,  because  a small  forcing  function  can  ex- 
plain the  amplitudes  in  the  stratosphere  and  tropo- 
sphere. However,  an  instability  consistent  with  the 
present  CISK  parameterization  should  eventually 
amplify  the  waves  to  a finite  amplitude  ih  comparable 
to  that  observed  for  large-scale,  cumulus-associated 
tropical  motions,  but  with  a temperature  amplitude 
too  large  to  be  realistic.  The  small  low-level  vertical 
velocities  of  these  waves  is  also  an  indication  that 
CISK  is  unlikely  to  be  important.  On  the  other  hand, 
the  tropospheric  amplitude  of  the  synoptic-scale  waves 
may  quite  possibly  be  supported  by  a CISK-type 
process.  These  waves  have  a wavelength  ~200(>-4000 
km  and  a period  ~5  days  which  lead  to  a phase  speed 
relative  to  the  ground  of  ~5-9  m s”‘.  Such  a phase 


’Lindzen  el  al.  (1975)  have  alternatively  considered  that  the 
CLSK  is  not  a lifting  process  but  a low-level  organizing  process 
which  does  not  require  lifting. 


speed  would  make  the  waves  greatly  attenuated  before 
reaching  the  stratosphere,  because  the  Doppler-shifted 
phase  speed  is  very  small  near  and  below  the  tropo- 
pause  and  the  damping  can  become  very  efficient 
(Lindzen,  1971). 
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ABSTRACT 

The  proljlem  of  scale-selection  of  Kelvin  waves  in  the  stratosphere  by  forcing  from  tropospheric  heating 
is  analyzed  using  a simple  linear  model.  The  effect  of  vertical  wind  shear  is  e.vciuded  because  the  phase 
speed  of  the  waves  is  fast  relative  to  the  range  of  the  mean  zonal  wind  in  the  vicinity  of  the  tropopause  at 
which  level  the  upward  energy  flu.v  due  to  forcing  is  evaluated.  Results  of  this  analysis  modify  Holton’s 
( 1973)  theory  in  that  1 ) the  forcing  is  most  efficient  for  the  longest  zonal  wavelength  even  if  the  heat  sources 
are  distributed  randomlv,  and  2)  the  most  favored  vertical  wavelength  of  the  e.’ccited  waves  is  about  twice 
the  vertical  scale  of  heating.  The  calculated  vertical  wavelengths  exceed  slightly  those  observed  and  the 
iliscrepancies  are  discussed. 


1.  Introduction 

There  is  ample  evidence  on  the  existence  of  atmo- 
spheric Kelvin  waves  in  the  etjuatorial  stratosphere 
(Wallace  ami  Koiisky.  l'>68;  and  others).  These  waves 
propagate  eastward  .ind  downward  with  a zonal  wave- 
number  of  1-2,  a periodii  itx  of  in- l.S  da\s,  and  a 
Vertical  wavelength  of  8-12  km  Their  zonal  wind 
component  has  a typical  amplitude  of  ID  m s~‘  at  the 
equator  while  no  appreciable  lluctuation  of  the  me- 
ridional wind  component  is  observed.  The  zonal  wind 
and  pressure  perturbations  are  both  symmetric  with 
respect  to  the  equator  and  are  close  to  geostrophic 
balance  in  the  meridional  direction,  but  the  waves 
behave  like  internal  gravity  waves  in  the  zonal  and 
vertical  directions.  .More  detailed  observational  and 
theoretical  descriptions  of  these  waves  are  given  in 
Wallace  (19731  and  Holton  (19751,  These  waves  are 
believed  to  play  a very  important  role  in  the  quasi- 
biennial oscillation  of  the  equatorial  stratosphere  by 
supplying  westerly  momentum  to  the  zonal  mean  flow 
(Lindzen  and  Holton,  1968;  Holton  and  Lindzen, 
1972). 

The  vertical  structure  of  the  Kelvin  waves  indicates 
that  they  are  maintained  by  energy  propagating  from 
the  troposphere.  In  addition,  the  possibility  that  they 
are  forced  by  interaction  with  mid-latitude  motions, 
as  suggested  by  -Mak  (1969),  is  quite  small  because 
of  the  small  meridional  wind  perturbations  which  rule 
out  the  importance  of  lateral  energy  flu.xes.  Dynamic 
instability  also  appears  unlikely  due  to  the  very  long 
zonal  scale  and  the  fast  westerly  phase  speed  of  the 
waves. 

Due  to  the  presence  of  the  large  amount  of  deep 
cumulus  omvection  in  the  tropical  troposphere,  it  has 


been  suggested  that  a jtlausible  energy  source  for 
these  waves  is  the  latent  heat  released  by  cumulus 
towers.  There  are  two  possible  mechanisms  by  which 
the  latent  heating  in  the  troposphere  may  maintain 
the  Kelvin  waves: 

1 ) two-way  interaction  such  that  the  cumulus 
convection,  which  supplies  energy  to  the  waves,  is 
itself  controlled,  by  the  large-scale  wave  motion  field. 
This  scale  interaction  results  in  an  unstable  situation 
which  is  generally  called  the  conditional  instability  of 
the  second  kind  (CISK  ' 

2)  The  waves  are  forced  by  cumulus  heating  which 
is  controlled  by  processes  unrelated  to  the  waves. 
Hayashi  (1970)  and  Lindzen  (1974)  have  investigated 
the  (irst  mechanism  and  found  that  the  growth  rates 
of  CISK  are  unsatisfactory  in  explaining  the  Kelvin 
waves.  Holton  (1972,  1973)  studied  the  second  mecha- 
nism with  a numerical  diagnostic  model  in  which  the 
diabatic  heating  due  to  cumulus  convection  is  speci- 
fied. Damping  in  the  form  of  Rayleigh  friction  and 
Newtonian  cooling  are  also  included.  He  found  that 
the  wave  response  in  the  stratosphere  to  a tropo- 
spheric heat  source  resembling  the  Kelvin  wave  mode 
has  a structure  in  close  agreement  with  observations. 
However,  the  observed  periodicity  of  the  Kelvin 
waves  must  be  specified  although  no  corresponding 
peak  in  the  cloud  brightness  spectra  can  be  found. 
On  the  other  hand,  his  numerical  calculations  suggest 
that  the  vertical  scale  of  the  excited  waves  is  close 
to  that  of  the  specified  forcing.  In  equatorial  wave 
theory  (Holton  and  Lindzen,  19681  the  vertical  scale 
determines  the  wavenumber-frequency  relationship. 
By  postulating  a red-noise  spectral  distribution  of  the 
tropospheric  heat  sources,  Holton  suggested  that  the 
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vertical  scale  and  the  large  low-frequency  variance 
of  the  forcing  result  in  a band-pass  selectivity  of  the 
Kelvin  waves. 

Holton’s  theory  may  indeed  be  valid,  especially 
because  an  oscillation  of  the  monsoon  heating  in  South 
Asia,  with  a quasi-periodicity  of  about  15  days,  has 
been  observed  during  the  northern  summer,  although 
it  is  not  clear  whether  this  periodicity  e.xists  in  other 
seasons.  The  purpose  of  this  paper  is  to  point  out 
that  the  atmosphere  acts  like  a band-pass  filter  even 
if  the  heat  sources  of  the  troposphere  are  distributed 
randomly  in  the  frequency  domain.  The  relationship 
between  the  vertical  scale  of  the  response  and  that 
of  heating  will  also  be  e.xamined. 

2.  The  forced  model 

The  zonal  mean  wind  in  the  tropical  stratosphere 
has  considerable  vertical  shear  which  may  influence 
the  upward  propagation  of  wave  energy  (Lindzen, 
1971).  However,  due  to  the  fast  phase  speed  of  the 
Kelvin  waves  the  shear  is  probably  not  very  signifi- 
cant at  the  tropopause  level  for  which  our  discussion 
will  be  most  relevant.  The  zonal  mean  wind  will 
therefore  be  neglected  in  our  simple  analysis.  The 
effect  of  damping  will  also  be  e.xcluded.  The  linearized 
zonal  momentum,  meridional  momentum,  hydrostatic, 
thermodynamic  energy  and  continuity  equations  on 
an  equatorial  beta-plane,  with  a single  zonal  wave- 
number  k,  a Doppler-shifted  phase  speed  c,  and  no 
perturbation  in  the  meridional  velocity  may  then  be 
written 


— ikcn=—ik<t>,  (1) 

^yU= , (2) 

dy 

d<t>  RT 

(3) 

dz  H 
Q 

— ikcT+U'T = — , (4) 

Cp 

d 


iku-\-e‘'" — (e~''"te)=0,  (5) 

dz 

where  «,  w,  T,  0 and  Q are  the  perturbation  zonal 
velocity,  vertical  velocity,  temperature,  geopotential 
and  the  diabatic  heating  rate,  respectively;  H is  a. 
constant  scale  height,  T the  static  stability,  Cp  the 
specific  heat  at  constant  pressure,  R the  gas  constant, 
li  the  meridional  gradient  of  the  vertical  component 
of  earth’s  vorticity,  y the  meridional  coordinate,  and 
z=—H  ln(/»  pn)  is  the  vertical  coordinate  with  p the 
pressure  and  pn  a reference  pressure. 

Eqs.  (l)-(2)  give  the  meridional  structure  of  the 


Kelvin  waves; 


«,  Ip  ® exp 


(6) 


where  the  condition  c>0  is  required  to  satisfy  the 
trapping  condition  on  an  equatorial  /J-plane.  If  a 
height-dependence  factor,  e\p{z/2H),  is  separated  from 
the  perturbation  quantities,  Eqs.  (1)  and  (3)-(5)  may 
be  combined  to  form  a single  equation  in  w : 


where 


ffhe'  O' 

i-av=— , 

as2  c2 

5 1 

4H'- 

R 

0'= 

CpII 

R 

5=-r 

// 


(7) 


(8) 


The  parameter  X is  a measure  of  the  vertical  wave- 
number.  If  the  heating  function  Q'  is  given,  (7)  can 
be  solved  with  suitable  boundary  conditions;  In  our 
problem  the  following  boundary  conditions  are  used : 


(0,  at  z=0  (9a) 

w=  I 

at  s=Z(  (9b) 

where 

Ci  = rC2. 


Here  s,  is  the  height  of  tropopause  and  the  condi- 
tion (9b)  results  from  the  requirement  that  latent 
heating  vanishes  at  z,.  The  parameter  r is  a reflection 
coefficient  which,  in  the  absence  of  vertical  wind  shear, 
is  given  by  the  model  specification  of  the  static  sta- 
bility distribution.  The  solution  to  (7)  at  z = 3i  may 
now  be  written  as 


w{z,)  = / sin(Xz)— rfz,  (10) 

X(l+r)  c‘ 

where  z^  is  the  height  of  the  cloud  base. 

3.  The  energy  flux  at  tropopause 

We  will  consider  the  vertical  energy  flu.x  0'a'', 
where  the  overbar  denotes  the  zonal  average  and 
0'=0e.\p(  — z 2/7),  at  the  tropopause  level  (z,),  to  be 
a measure  of  the  etliciency  of  forcing.  Despite  the 
actual  existence  of  the  vertical  shear  of  the  zonal 
mean  flow  in  the  stratosphere,  the  response  of  the 
waves  should  be  largely  determined  by  this  quantity. 
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At  the  solution  is  given  by  (9b)  and  only  the  Cj 
term  needs  to  be  considered  for  upward  energy  propa- 
gation or  downward  phase  propagation. 

Eqs.  (1)  and  (5)  may  now  be  rewritten  as 

— ikcu'~—ik<i>', 
iku'—^iX-l ^w'=0, 

where  m'  = m exp(— z/2/7).  Front  these  relationships  it 
follows  immediately  that 


Since  //»7  km  for  the  atmosphere,  /\»1  (2//)  if  the 
vertical  wavelength  «88  km.  This  is  the  case  for  all 
practical  purposes  so  we  may  write 

cA 

</>'==— te'. 

k 

From  the  detinilion  of  .\  given  in  Eq.  (8)  we  may  also 
write 

c\^S\ 

which  leads  to 

51 


and 

51_ 

(11) 

k 

Eq.  (11)  indicates  that,  if  the  amplitude  of  the  ver- 
tical velocity  perturbation  at  the  tropopause  level  is 
ti.xed,  the  upward  energy  fiu.x  will  be  inversely  propor- 
tional to  the  zonal  wavenumber.  E.xamination  of  (7) 
reveals  that  w'  depends  only  on  X (c  is  uniquely  de- 
termined by  X)  as  long  as  Q is  independent  of  k.  Thus 
for  given  vertical  scale  the  forcing  at  the  tropopause 
level  becomes  maximum  for  the  longest  horizontal 
wavelength.  This  result  obviously  remains  the  same 
if  one  considers  the  upward  momentum  flux  ti'w', 
because  u'  = 4ifc. 

In  order  to  e.xamine  the  selection  of  vertical  scales, 
we  assume  a heating  function 

|m5(/(y)sinir( ),  z'^z^Zc 

Q'  = l \ Az  / (12) 

lo,  z>z,  or  z<z» 

where  m specifies  the  strength  of  the  heating  and 
Ac  = z,— z,  is  a measure  of  the  vertical  scale  of  the 
heating.  The  function  q(y)  which  specifies  the  me- 
ridional distribution  of  heating  will  be  assumed  to 


take  the  form  of  (6)  for  the  time  being.  The  vertical 
profile  as  given  by  (12)  resembles  closely  that  of  the 
typical  observations  (Chang,  1976).  For  the  simplest 
case,  we  shall  first  assume  that  the  static  stability 
parameter  5 has  a constant  tropospheric  value 

5 = 5, 

throughout  the  atmosphere.  In  this  case  r=0  and  only 
upward  energy  propagation  at  z=z,  is  allowed.  The 
solution  (10)  at  z,  is  now  found  to  be 

w'=wexp[— iXz,— /3y'/(2c)],  z = z„  (13) 

where 


mwAz  X'-f  1/(477)2 

w= 

X ir2-X2Az2 


sin  -(z,+Zr) 
12 


(14) 


denotes  the  amplitude  of  the  solution  in  this  simple 
model.  We  note  here  that  no  actual  singularity  exists 
in  (14).  The  vertical  energy  flux  at  the  equator  and 
z=Z(  is  thus 

5‘_ 

0'w'= — w~. 

2k 


From  (14)  if  one  assumes  X»1  (477^),  it  can  be  shown 
that  this  energy  flux  reaches  maximum  at 


where 


TT 

X=— , 
Az 


mV  r V 

(Fa sin'd (s,-)-z,)  . 

16  L2Az  J 


The  tropospheric  forcing  is  therefore  most  erficient  at 
a vertical  wavelength 

2ir 

L,=  — = 2Az, 

X 


or  twice  the  vertical  scale  of  the  heating.  This  result 
is  somewhat  similar  to  those  of  Green  (1965)  and 
Lindzen  (1966).  A plot  of  the  energy  flux  profile  at  z, 
is  given  in  Fig.  1 in  which  it  is  assumed  that  z,  = 15  z^ 
and  H = lZc.  If  we  measure  the  vertical  scale  of  Hol- 
ton’s (1972,  1973)  heating  function  the  value  Az  can 
be  found  to  be  ~11  km,  which  predicts  a maximum 
response  of  22  km  vertical  wavelength  in  the  tropo- 
sphere. Since  in  reality  the  stratospheric  static  sta- 
bility is  about  three  times  that  of  the  troposphere, 
using  (8)  we  may  estimate  that  the  stratospheric 
vertical  wavelength  corresponding  to  a 22  km  tropo- 
spheric value  would  be 


= 22X3-‘«12.8  km. 


where  the  subscript  s denotes  stratosphere.  If  we 
assume  Az«  14  km,  then  7„«  16.3  km. 
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It  may  be  of  some  interest  to  consider  the  effect 
of  different  static  stability  between  stratosphere  and 
troposphere  on  the  solution  (13),  because  the  “effi- 
ciency” of  reflection  of  wave  energy  at  2 = 3,  may  be 
scale-dependent.  In  this  case  we  have 

S=S„  for  2>2(, 

S=St,  for  2^  2,. 

.\nother  boundary  condition  is  now  needed  and  the 
simplest  form  available  is  the  radiation  condition 
above  z,,  i.e., 

U''  = A exp[— *X,2— /3)r/(2c)],  for  z>z,,  (15) 

where  X,  is  the  vertical  wavenumber  in  the  strato- 
sphere. The  continuity  conditions  of  w'  and  <t>'  across 
2(  lead  to 

l-(5./5.)‘ 

r = exp(  — 2/X,2,),  for  S.>S,, 

l + (5./5,)‘ 

where  X,  is  the  vertical  wavenumber  in  the  tropo- 
sphere. The  solution  (10)  at  2,  in  the  form  of  (15)  is 

2X(  exp[t(X,— Xi)2jil' 

A , 

(X.-fX,)  — (X.— X,)  exp(— 2jXt2,) 

where  w is  given  by  (14)  with  all  X replaced  by  X,. 
The  vertical  energy  flux  at  the  equator  and  2,  is  now 

5‘ 

(i>V^—AA*==>—Ew\ 


E=2l{l+S./Sd+{i-S./S,)  cos2X,2,]-‘ 

is  the  effect  of  different  static  stabilities  above  and  below 
2i  and  the  asterisk  denotes  complex  conjugate.  For 
S,=  3Si,  E varies  between  a maximum  of  1 when 
X(=0,  r/zt,  2ir/zi,  ....  and  a minimum  of  0.33  when 

Xi=t/(22(),  .3t/(22(),  St/(23(),  Assuming  2(=  15 

km  and  A2=  14  km,  the  variation  of  E is  plotted  in 
the  upper  section  of  Fig.  1 and  the  modified  energy 
flux  profile  is  indicated  by  the  dashed  curve  in  the 
lower  section.  It  is  seen  that  the  most  efficient  forcing 
is  slightly  larger  than  twice  the  scale  of  heating.  Hence 
we  may  conclude  that  the  variation  in  static  stability 
between  troposphere  and  stratosphere  does  not  cause 
great  effect  in  the  vertical  scale  of  the  most  favored 
response. 

4.  Concluding  remarks 

We  have  shown  that  for  randomly  distributed 
tropospheric  heat  sources  the  Kelvin  waves  of  the 


Li- 


Fig.  1.  Profile  of  vertiLal  wave  energy  flux  as  a function  of 
vertical  wavelength  L,.  The  solid  curve  in  the  lower  diagram  is 
for  the  case  5, =5,  and  the  dashed  curve  for  the  case  5,-35,.  The 
upper  diagram  is  the  effect  on  E of  the  difference  in  static  stability 
on  the  vertical  wave  energy  flux  profile  for  the  case  5,—35i. 

largest  zonal  scale  are  excited  most  efficiently  in  the 
absence  of  damping.  This  may  explain  why  most  of 
the  Kelvin  waves  observed  in  the  lower  stratosphere 
are  of  zonal  wavenumber  1 and  some  are  of  wave- 
number  2.  The  vertical  wind  shear  has  been  excluded 
from  our  analysis  but  it  should  not  significantly  alter 
this  conclusion. 

The  finding  that  the  most  favored  vertical  wave- 
length is  about  twice  the  vertical  scale  of  heating 
gives  a reasonable  estimate  for  Kelvin  waves  in  the 
lower  stratosphere,  although  the  calculated  wave- 
length of  13-16  km  is  slightly  larger  than  those  ob- 
served. Several  influences  which  are  excluded  in  our 
model  may  alter  this  aspect.  For  example,  Lindzen 
(1971)  found  that  the  focusing  effect  of  the  vertical 
wind  shear  in  the  stratosphere  would  reduce  the 
vertical  scale  of  the  waves.  .\  two-scale  expansion 
technique  can  be  similarly  applied  to  our  forced  model 
and  the  resultant  prevailing  vertical  wavelength  in 
the  stratosphere  should  be  reduced. 

Another  effect  which  has  been  excluded  is  the 
actual  north-south  distribution  of  the  heating  func- 
tion. The  meridional  structure  equation  (6)  indicates 
that  the  meridional  scale  of  the  waves  is  coupled  with 
the  vertical  scale  through  the  phase  speed  c.  The 
variation  of  the  meridional  scale  with  the  vertical 
scale  may  be  considered  slow  because  the  e-folding 
width  of  the  Gaussian  distribution  is  proportional  only 
to  the  square  root  of  c,  which  is  approximately  pro- 
portional to  the  vertical  wavelength.  In  the  real 
atmosphere  the  meridional  heating  profile  q{y)  in  (12) 
does  not  adjust  to  different  vertical  scales  as  we  have 
assumed  in  our  calculations.  A proper  treatment  of 
this  problem  is,  perhaps,  to  multiply  the  previously 
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calculated  response  by  the  projection 

f q{y)  exp[-/JyV(2c)]<i’y 
J 0 

P= ^ . (16) 

f exp[-/Jyy(2c)]Vy 
J 0 

Here  we  could  specify  q{y)  as  we  did  for  the  vertical 
heating  profile,  although  there  appears  to  be  more 
observational  uncertainty.  For  the  purpose  of  the 
present  qualitative  discussion,  however,  this  is  not 
strictly  necessary.  We  know  that  the  latent  heating 
is  largely  confined  to  the  tropics,  so  that  the  meridional 
extent  of  q(y)  must  be  quite  limited.  Therefore,  the 
effect  of  (16)  should  be  such  that  waves  whose  vertical 
scales  are  larger  than  that  corresponding  to  the  me- 
ridional scale  of  q{y)  will  suffer  from  having  a large 
portion  of  their  north-south  domain  unsupported  by 
heating.  If  we  assume  that  the  meridional  scale  of  q{y) 
can  be  represented  by  an  e-folding  width  yd«  1500  km 
of  a Gaussian  distribution  about  the  equator,  and 
that  5(=  1.22X10“*  5“-  which  corresponds  to  r = 3 K 
km“‘,  then  those  waves  with  a tropospheric  vertical 
wavelength 

L,>ir;iSt~^yd-~  l-i-2  km, 
or  a stratospheric  wavelength 

Z.,>8.2  km, 

will  have  their  amplitude  reduced  by  (16).  The  reduc- 
tion will  be  more  severe  for  longer  vertical  wave- 
lengths. This  obviously  would  also  shorten  the  pre- 
viously calculated  most  favored  vertical  scales  of  the 
Kelvin  waves. 
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In  a recent  paper  by  Kuo  (1975),  a condition  was 
used  in  a CISK  parameterization  to  e.xplain  the  time 
and  length  scales  of  the  observed  mixed  Rossby -gravity 


waves  and  the  Kelvin  waves  in  the  tropics.  This  condi- 
tion, as  stated  by  Kuo  on  p.  2231,  is  that  “the  depletion 
of  moisture  through  precipitation  must  be  replenished 
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by  evaporation  from  the  surface  during  a complete 
period  of  the  disturbance.”  Eq.  (15b)  in  his  paper 
specifies  this  condition  as  follows : 

r (£+rate  of  moisture  convergence)  ^ maximum  total 
precipitation  given  by  deep  cumulus  convection 
during  a period, 

where  t is  the  period  of  the  wave  disturbances  and  E 
the  evaporation  rate.  On  p.  2234  he  derived  the  limiting 
period  to  for  which  the  equal  sign  in  (15b)  holds.  This 
is  given  by  his  Eq.  (23) : 


P, 


where  Pt  is  the  average  maximum  total  precipitation 
over  a period.  This  maximum  total  precipitation, 
according  to  the  argument  on  p.  2232,  is  always  pro- 
duced by  the  disturbance  “once  the  convective  insta- 
bility is  activated  by  the  large-scale  flow,”  because  the 
“convection  will  proceed  as  an  overturning  process  and 
produce  the  maximum  amount  of  precipitation  possible 
for  the  given  air  mass.”  The  value  of  Pt  was  estimated 
by  subtracting  the  vapor  contents  of  the  outflow  column 
from  the  inflow  column  as  was  schematically  illustrated 
in  his  Fig.  3.  The  computed  tq  was  then  used  as  a cutoff 
period  in  the  CISK  parameterization  [Eq.  (24)]  in  such 
a way  that  the  heating  function  is  proportional  to  the 
factor  r/ro  if  r<To.  The  implication  of  his  analysis  is 
that,  “disturbances  with  periods  <to  cannot  be  sup- 
ported fully  by  deep  cumulus  convection  initiated  by 
such  disturbances  because  the  moisture  content  of  the 
atmosphere  will  be  depleted  by  the  convective  process.” 

The  above  arguments  are,  however,  inconsistent  with 
an  essential  part  of  the  CISK  parameterization  used  in 
the  paper,  namely  that  the  large-scale  heating  rate  is  pro- 
portional to  the  vertical  velocity  IFj  at  a low  reference 
level  Z4  when  IFi>0.  In  this  parameterization  when  the 
wave  amplitude  at  an  early  stage  of  the  development  is 
small,  W t,  is  small  and  so  is  the  heating  rate,  which 
implies  that  the  precipitation  rate  is  small.  As  the 
waves  grow,  Wi,  increases  and  the  precipitation  rate  in 
turn  increases.  Since  the  total  precipitation  of  a 


continuous  wave  chain  cannot  e.xceed  the  total  evapora- 
tion in  the  tropical  belt,  the  restraining  effect  is  such 
that  the  precipitation  rale  must  not  exceed  the  evapora- 
tion rate,  or  3.75  g c m"’  day*  according  to  the  calcula- 
tions quoted  bv  Kuo.  Therefore  the  total  moisture 
availability  in  the  tropics  should  constrain  the  maxi- 
mum amplitude  the  waves  can  grow  to,  rather  than  the 
period  of  the  waves.  The  intensity  of  the  cumulus-scale 
convection  is  determined  by  the  large-scale  moisture 
convergence  below  Zi,  in  the  CISK  theory  because  this 
convergence  provides  the  high  8e  for  the  ascending  air. 
For  more  sophisticated  parameterization  theories,  such 
as  that  of  Arakawa  and  Schubert  (1974)  which  is  a 
form  of  the  convective  adjustment,  the  same  incon- 
sistency would  occur.  In  these  theories  the  effect  of 
cumulus  convection  is  also  entirely  a response  to  the 
large-scale  forcing  and  therefore  its  intensity,  as  re- 
flected by  the  net  large-scale  heating,  should  be  con- 
trolled by  the  large-scale  flow. 

On  the  other  hand,  another  condition  mentioned  by 
Kuo  [Eq.  (15a)]  is  a valid  one,  and  may  be  quite 
important  in  selecting  the  wave  scales.  This  condition 
requires  “a  disturbance  of  sufficiently  long  duration  and 
high  intensity  able  to  lift  the  lower  atmosphere  above 
the  condensation  level,  thereby  activating  the  deep 
cumulus  convection.”  In  other  words,  waves  with  short 
periods  are  unlikely  to  grow.  Kuo’s  Eq.  (24),  which 
contains  a cutoff  period  may,  in  a way,  be  viewed  as  a 
crude  representation  of  this  condition  although  the 
precise  requirement  can  only  be  treated  as  a nonlinear 
problem.  Based  on  this  condition  alone,  waves  whose 
Doppler-shifted  phase  speed  equals  zero  for  are 
most  favored  for  development.  This  seems  to  be  the 
case  in  the  tropical  western  Pacific  where  the  usually 
observed  warm-core  synoptic -scale  waves  have  a 
typical  easterly  phase  speed  of  5-8  m s“‘,  quite  close  to 
the  low-level  mean  zonal  wind. 
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ABSTRACT 

In  this  paper  we  deal  with  the  interpretation  of  observed  oscillations  in  the  tropical  troposphere  and 
stratosphere  within  the  framework  of  the  equatorial  wave  theory.  A difficulty  with  this  problem  arises 
when  one  compares  the  short  vertical  wavelength  (or  equivalent  depth)  predicted  by  the  classical  theory 
and  the  observed  large  vertical  scales  associated  with  the  low  Uoppler-shifted  frequencies  of  the  tropospheric 
oscillations.  In  this  analysis  it  is  shown  that  the  inclusion  of  simple  linear  damping,  justified  by  budget  studies 
which  revealed  the  important  role  of  cumulus  momentum  transport,  has  a strong  influence  at  low  fre- 
quencies on  the  forced  equatorial  waves  and  results  in  two  types  of  dispersive  relationships.  The  first  type 
is  characteristic  of  the  regular  internal  gravity  waves  which  have  fast  phase  speeds  and  weak  vertical 
attenuation.  The  second  type  is  dominated  by  the  viscous  damping  time  scale  and  has  slow  phase  speeds 
and  strong  vertical  trapping.  The  theory  predicts  that  the  stratospheric  oscillations  may  be  identified  with 
the  first  type  and  the  tropospheric  oscillations  with  the  second.  In  the  case  of  Kelvin  waves  the  results 
can  be  used  to  explain  consistently  both  the  observed  stratospheric  Kelvin  waves  and  the  planetary-scale 
Kelvin-like  oscillations  in  the  troposphere  including  the  40-50  day  oscillation  and  the  monsoon  and  Walker 
circulations.  Possible  implications  with  respect  to  other  waves  in  the  tropics  are  also  discussed. 


1.  Introduction 

Based  on  spectral  and  cross-sp)ectral  analyses  of 
long  time  series  (5-10  years)  of  tropical  station  data, 
Madden  and  Julian  (1971,  1972)  detected  a global- 
scale  oscillation  in  the  tropical  trofiosphere  with  a 
periodicity  of  40-50  days.  A degree  of  stationarity 
of  this  oscillation  in  time  was  also  established  by 
them  by  examining  the  limited  amount  of  station 
pressure  data  available  from  the  1890’s.  Their  results 
indicate  that  this  oscillation  is  most  prominent  in  the 
zonal  velocity  component  and  the  surface  pressure, 
with  the  maximum  amplitude  situated  along  the 
equator.  There  are  some  time  and  spatial  variations 
within  a 40-50  day  cycle  as  depicted  in  Fig.  1 which 
is  reproduced  from  a schematic  diagram  constructed 
by  Madden  and  Julian  (1972),  but  basically  the  oscil- 
lation propagates  eastward  at  a fairly  steady  phase 
speed  and,  at  most  times,  exhibits  a zonal  wave- 
number  1 structure.  Madden  and  Julian  (1972)  argued 
that  enhanced  deep  cumulus  convection  is  usually 
associated  with  large-scale  upward  motion  and  the 
oscillation  appears  as  two  vertical  circulation  cells 
along  the  equatorial  zonal  plane. 

The  above  structure  resembles  that  of  a theoretical 
atmospheric  Kelvin  wave  (Holton  and  Lindzen,  1968) 
in  several  ways:  the  direction  of  propagation,  the 
meridional  structure  and  the  absence  of  large  amplitude 
in  the  meridional  wind  component  However,  for  a 


given  phase  speed  as  slow  as  that  observed,  Kelvin 
wave  theory  would  predict  a fairly  short  vertical 
wavelength  which  does  not  appear  to  agree  with  that 
of  the  observed  oscillation.  This  discrepancy,  which 
has  been  a major  difficulty  for  the  theoretical  under- 
standing of  the  40-50  day  oscillations,  has  also  been 
encountered  by  several  meteorologists  in  their  efforts 
to  interpret  various  types  of  large-scale  tropical  wave 
motions  using  the  equatorial  wave  theory.  In  a linear 
analysis  Lindzen  (1967)  has  shown  that  for  synoptic 
(and  larger)  time  scales,  waves  on  an  equatorial 
beta-plane  generally  propagate  in  the  vertical  as  in- 
ternal gravity  waves  with  rather  short  vertical  wave- 
lengths. This  result  was  evident  in  a study  of  the 
response  of  the  tropical  atmosphere  to  stationary 
heating  by  Webster  (1972,  1973).  In  a two-level 
(750  and  250  mb)  numerical  model  he  found  that  the 
equatorial  response  in  an  easterly  basic  current  has 
a horizontal  structure  that  is  characteristic  of  Kelvin 
waves.  However,  his  analytical  solution  for  the  case 
of  an  idealized  basic  flow  possesses  very  short  vertical 
wavelength  [~0(1  km)]  and  appears  to  suggest  that 
higher  vertical  resolution  is  required  for  the  nu- 
merical model. 

The  short  vertical  wavelength  predicted  by  the 
equatorial  wave  theory  arises  from  the  internal  gravity 
wave  character  of  these  waves.  Since  for  pure  gravity 
waves  the  phase  speed  c is  given  by  c=(gA)*,  where 
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Ftc.  1.  Schematic  depiction  of  the  time  and  space  (zonal 
plane)  variations  of  the  circulation  celb  associated  with  the 
40-50  day  oscillation.  The  mean  pressure  disturbance  is  plotted 
at  the  bottom  of  each  chart  with  negative  anomalies  shaded. 
Regions  of  enhanced  large-scale  convection  are  indicated  sche- 
matically by  the  cumulus  and  cumulonimbus  clouds.  The  relative 
tropopause  height  is  indicated  at  the  top  of  each  chart  (From 
Madden  and  Julian,  1972.) 

g is  the  gravitational  constant  and  k a depth  scale 
for  the  motion  (usually  called  the  equivalent  depth), 
the  vertical  wavelength  tends  to  Le  proportional  to 


the  phase  speed.  Within  the  range  of  large-scale 
tropical  motions  in  the  tropical  troposphere  the  phase 
speed  is  usually  slow  compared  to  the  external  gravity 
wave  phase  speed,  so  the  vertical  wavelength  is  con- 
sequentially small  compared  to  that  given  by  the 
scale  height.  The  observed  Kelvin  waves  (Wallace 
and  Kousky,  1968)  and  mixed  Rossby-gravity  waves 
(Yanai  and  Maruyama,  1966)  in  the  lower  strato- 
sphere have  rather  large  Doppler-shifted  phase  speeds 
(relative  to  the  mean  zonal  wind)  so  that  their  ob- 
served vertical  wavelengths  are  in  better  agreement 
with  the  theory.  In  Table  1 we  summarize  the  com- 
parison between  the  observed  vertical  wavelengths 
and  that  given  by  the  equatorial  wave  theorj'  for 
three  major  zonal  wind  oscillations  in  the  tropics: 
the  stratospheric  Kelvin  wave,  the  40-50  day  oscilla- 
tion and  the  planetary-scale  stationary  motion  which 
could  include  the  east-west  overturning  Walker  cir- 
culation along  the  equatorial  plane  (Bjerknes,  1969; 
Webster,  1973)  or  the  winter  monsoon  (Krishnamurti 
el  al.,  1973;  Webster,  1973).  The  theoretical  vertical 
wavelength  is  calculated  using  the  Doppler-shifted 
phase  speed  with  the  mean  zonal  flow  specified  by 
an  averaged  tropospheric  easterly  speed  of  — 5 m s~‘. 
This  specification  is  based  on  the  evidence  that  all 
these  oscillations,  including  the  stratospheric  Kelvin 
wave,  are  generated  by  energy  sources  in  the 
trojxisphere. 

The  equatorial  wave  theory  used  in  Table  1 is 
basically  inviscid  because  it  does  not  include  a sig- 
nificant frictional  effect  in  the  free  atmosphere.  In  a 
numerical  diagnostic  study  of  the  seasonal-mean  cir- 
culation at  200  mb  during  the  northern  summer, 
Holton  and  Colton  (1972)  found  that  a very  strong 
damping  term  (corresponding  to  linear  damping  time 
scale  of  ~ 1 day)  is  needed  in  the  linearized  barotropic 
vorticity  equation  to  balance  the  generation  of  plane- 
tary-scale vorticity  by  horizontal  divergence.  Since 
the  planetary-scale  divergence  field  is  largely  a result 
of  the  strong  condensation  heating  in  the  middle 
troposphere  due  to  the  vigorous  cumulus  convection 
associated  with  the  summer  monsoon,  Holton  and 
Colton  hypothesized  that  the  vertical  mixing  of 
vorticity  and  momentum  by  cumulus  convection 
could  account  for  this  strong  damping.  This  hypothesis 
is  also  consistent  with  several  budget  studies  of  syn- 
optic-scale, cumulus-related  tropical  motions  (Reed 
and  Recker,  1971;  Williams  and  Gray,  1973;  Reed 
and  Johnson,  1974;  Chu,  1976;  and  others),  and  the 
numerical  diagnostic  study  by  Chang  el  al.  (1975). 
Recently  Stevens  el  al.  (1977)  also  show  that  this 
strong  damping  mechanism  is  responsible  for  reducing 
the  theoretically  calculated  amplitude  of  the  tropical 
waves  forced  by  heating  to  the  observed  magnitude. 
Thus  the  inclusion  of  this  mechanism  in  the  equatorial 
wave  theory  seems  to  be  necessary  for  the  proper 
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T ABL£  1 . Comparison  of  observed  tropical  zonal  oscillations  and  the  inviscid  equatorial  wave  theory.  The  theoretical  vertical  wavelengths 
are  based  on  the  observed  Doppler-shifted  phase  speeds  which  assume  a basic  zonal  wind  of  d»  — 5 m s~‘.  See  text  for  details. 


Observations 
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Oscillations 

Period 
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phase  speed 
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(m  s-*) 
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(km) 

Stratospheric  Kelvin  waves 
(Wallace  and  Kousky, 

10-15 

1-2 
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1968) 

40-50  day  oscillations 
(Madden  and  Julian, 

40-50 

1 

9-11 

14-16 

15-30 

8-9 

1972) 

Stationary  pianetary 
circulations  (monsoon 
and  Walker) 

00 

1-2 

0 

5 

15-30 

3 

• u » — 5 m s”'. 


interpretation  of  the  observed  tropical  motions  forced 
by  cumulus  heating. 

The  purpose  of  this  paper  is  to  study  the  effect 
of  the  strong  damping  mechanism  called  for  by  pre- 
vious budget  studies  on  tropical  wave  dynamics.  The 
simplest  form  of  linear  drag  will  be  used  to  represent 
this  effect  as  was  done  by  Holton  and  Colton  (1972). 
We  will  show  that  this  effect  in  some  cases  is  more 
than  just  an  expected  reduction  in  the  amplitude  of 
the  thermally  forced  waves  or  a near-trivial  correction 
in  the  frequency  equation.  Moreover,  in  these  cases 
it  can  account  for  most  of  the  discrepancies  in  Table  1 
and  provide  a consistent  theoretical  basis  for  the 
interpretation  of  the  tropical  oscillations  in  both  the 
stratosphere  and  the  trofwsphere. 

2.  Basic  equations 

Since  the  vertical  structure  equation  of  all  equatorial 
waves  is  identical  to  that  for  two-dimensional  gravity 
waves,  we  will  consider  only  the  Kelvin  wave  case 
for  the  purpose  of  simplicity.  Setting  the  meridional 
velocity  zero  we  may  write  the  linearized  zonal  mo- 
mentum, meridional  momentum,  hydrostatic,  thermo- 
dynamic energy  and  continuity  equations  on  an 
equatorial  beta-plane,  with  a single  zonal  wave- 
number  k and  constant  phase  sp)eec  c relative  to  the 
basic  zonal  flow  as 

— ikcu  — —ik<^— Du,  ( 1 ) 

0yu= , (2) 

0y 

00  RT 


Q 

-ikcT+wT= DT,  (4) 

Cp 

0 

iku+e‘i»—(e-‘i»w)=0,  (5) 

dz 

where  «,  w,  T,  0 and  Q are  the  perturbation  zonal 
velocity,  vertical  velocity,  temperature,  geopotential 
and  the  diabatic  heating  rate,  respectively;  £7  is  a 
constant  scale  height,  F the  static  stability,  Cp  the 
specific  heat  at  constant  pressure,  R the  gas  constant, 
0 the  meridional  gradient  of  the  vertical  component 
of  earth’s  vorticity,  y the  meridional  coordinate  ; and 
z=  —H  \n(p/po)  is  the  vertical  coordinate  with  p the 
pressure  and  po  a reference  pressure.  The  damping 
effect  is  represented  by  the  linear  drag  coefficient  D in 
the  zonal  momentum  equation.  A Newtonian  cooling 
term  with  the  same  coefficient  is  included  in  the 
thermodynamic  energy  equeUon  in  order  to  facilitate 
the  analysis. 

Eqs.  (l)-(2)  give  the  meridional  structure  of  the 
Kelvin  waves ; 

M,0<xexp(-;8y*/22),  (6) 

where 

and  the  condition  Re(^)  = c>0  is  required  to  satisfy 
the  trapping  condition  on  an  equatorial  /3-plane.  If  a 
height-dependence  factor,  exp(z/2H),  is  separated  from 
the  perturbation  quantities,  Eqs.  (1)  and  (3)-(5)  may 
be  combined  to  form  a single  equation  in  w : 


3b' 


Q' 

XV—, 

2* 


(7) 
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where 


R 

Q'= 

CfH 


(8) 


The  parameter  X is  a measure  of  the  vertical  wave- 
number.  If  the  heating  function  Q'  is  given,  (7)  can 
be  solved  with  suitable  boundary  conditions.  In  our 
problem  the  following  boundary  conditions  are  used: 


0,  at  3 = 0 (9a) 

w=  • 

at  3 = 3;  (9b) 

where 

Ci=rC2. 


Here  3,  is  the  height  of  tropopause  and  the  condi- 
tion (9b)  results  from  the  requirement  that  latent 
heating  vanishes  at  3,.  The  parameter  r is  a reflection 
coeflflcient  which,  in  the  absence  of  vertical  wind 
shear,  is  given  by  the  model  specification  of  the  static 
stability  distribution.  The  solution  to  (7)  may  be 
written  in  the  form 


■w{z)=- 


X(l-hr) 
sinX3  r“ 


—t\t  i*tt 

/ sinXs — dz,  3^3, 

r)  J,c  gh 


(10a) 


z r“ 

— I {re'^‘ 

r)  J. 

C [ sir 
J f« 


Q' 

— dz 

gh  X(l-hr) 


Q' 


X/  sinX3— ds,  3,>3>3c,  (10b) 

gh 


sinXs  /•“  Q' 

w{z)  = / — dz,  3c^3^0.  (10c) 

X(l-hr)  J gh 


Here  Zt  is  the  height  of  cloud  top  and  3^  the  height 
of  cloud  base.  Chang  (1976)  has  used  the  above 
equations  to  show  that,  in  the  case  of  Z7=0,  the 
vertical  energy  flux  of  Kelvin  waves  due  to  pressure 
work  is 

5‘_2 

(p'w'si Itl'  , (11) 

k 


and  Lindzen,  1966),  Chang  (1976)  was  able  to  show 
that  the  Kelvin  wave  response  given  by  the  inviscid 
theory  agrees  with  the  observed  stratospheric  Kelvin 
waves. 

If  we  consider  downward  phase  propagation  (up- 
ward energy  propagation)  only,  Re(X)>0.  For  H 
~7  km  and  vertical  wavelength  «88  kra,  Eq.  (8) 
may  be  approximated  by 

cX=5*.  (12) 

In  the  inviscid  case,  Z)=0  and  c=c  so  that  X is  pure 
real.  In  the  case  of  &=CT+iCi  where  Cr=r>0 

and  Ci=D/k>0;  thus  X becomes  complex,  i.e.,  X=Xr 
-fiXi.  Here  Xi  describes  the  exponential  increase 
(Xi>0)  or  decrease  (X,<0)  of  the  amplitude  with 
height.  Eq.  (12)  now  becomes  quadratic  and  has  the 
following  roots : 

5‘±(5-4X,V)‘ 

Cri  n = 5>0,  (13) 

2X, 

-5i±(5-4X,V)‘ 

X,.,„= <0.  (14) 

2c,' 


The  subscripts  I and  II  are  being  used  to  denote  the 
plus  and  minus  roots,  respectively.  It  is  obvious  that 
Cri>Crii  and  jX,,  | < j X,,[  | , hence  the  first  solution 
(which  shall  be  called  Mode  I)  has  a fast  horizontal 
phase  speed  and  can  propagate  more  efificiently  into 
higher  levels  while  the  second  solution  (Mode  II)  is 
slower  and  attenuates  faster  above  the  source  region. 
In  addition,  Eqs.  (13)  and  (14)  indicate  that  the  waves 
are  excitable  only  if 


Xr 


(15) 


so  that  the  effect  of  damping  causes  a short-wave 
cutoff  in  the  vertical  or  a long-wave  cutoff  in  the 
horizontal. 

If  we  take  the  inviscid  limit  of  D—>0,  Eqs.  (13) 
and  (14)  reduce  to 


c„=5VX,l 
X„=0  I’ 

Crii=0  I 
Xi'fi  = * ^ 


(16) 

(17) 


where  the  overbar  denotes  the  zonal  average  and 
(t>'  = <t>exp(~z/2H).  Eqs.  (11)  and  (7)  indicate  that 
for  this  case  the  upward  forcing  of  Kelvin  waves  is 
most  eificient  for  long  waves.  Augmented  by  the  fact 
that  waves  with  a vertical  scale  equal  to  twice  that 
of  forcing  are  most  efficiently  excited  (Green,  1965 


Comparing  (16)  and  (12)  it  is  seen  that  Mode  I cor- 
responds to  the  regular  inviscid  solution.  Substituting 
(17)  into  (l)-(5)  leads  to  a trivial  solution  of  vanish- 
ing wave  amplitude.  Thus  Mode  I may  be  called  the 
“regular  mode”  and  Mode  II  the  “viscous  mode,” 
because  the  latter  is  non-trivial  only  when  Dt^O. 
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Fig.  2.  Real  Doppler-shifted  phase  speed  c,  as  a function  of  vertical  wavelength  i,  and  damping 
time.  For  the  convenience  of  application  to  other  types  of  equatorial  waves,  the  abscissa  can  also  be 
designated  by  (g/i)r*i  where  A is  the  equivalent  depth. 


Fig.  2 shows  the  real  phase  speed  Cr  as  a function  3. 
of  vertical  wavelength  Z,(  = 2!r  \,)  and  the  viscous 
damping  time  scale  l/D.  For  D=0  (damping  time 
scale  « j the  dispersion  relationship  of  Mode  I is  a 
straight  line  of  Cr'^  L„  typical  of  the  internal  gravity 
waves.  In  addition,  the  vertical  coordinate  (c,=0) 
represents  the  trivial  solution  of  Mode  II.  For  Z)?sO 
the  two  modes  are  represented  by  the  quadratic 
curves.  The  effect  of  damping  on  Mode  I is  a rela- 
tively small  modification  with  L,  increasing  as  c,  in- 
creases. The  modification  for  Mode  II  is  more  sub- 
stantial and  L,  generally  decreases  as  c,  increases. 

The  two  modes  are  separated  by  the  line  c,  = Ci.  For 
Mode  I c,>c„  i.e.,  the  inertial  time  scale  is  faster 
than  the  damping  time  scale.  The  reverse  is  true  for 
Mode  II. 

The  imaginary  vertical  wavenumber  X,  is  shown  in 
Fig.  3 as  a function  of  vertical  wavelength  and  damping 
time  scale.  Here  the  magnitude  of  Xi  decreases  with 
increasing  L,  for  Mode  I and  increases  with  L,  for 
Mode  II.  The  c,=c,  curve  again  provides  the  parti- 
tion of  the  two  modes. 

From  Figs.  2 and  3 it  can  be  seen  that,  if  a certain 
vertical  wavelength  is  most  efficiently  excited  by  a 
given  forcing  function,  both  modes  will  be  excited 
for  Z?^0.  The  regular  mode  will  have  fast  zonal  phase 
sfjeed  and  less  vertical  attenuation,  while  the  viscous 
mode  will  have  slow  zonal  phase  speed  and  is  more 
trapped  in  the  vertical  away  from  the  source  region. 


Response  of  specified  heating 

The  Kelvin  wave  response  of  the  tropical  atmo- 
.sphere  can  be  calculated  if  the  vertical  heating  func- 
tion Q'  in  (7)  is  known.  As  in  Chang  (1976)  we  will 
use  the  vertical  heating  profile 

Q! <x Sq (y)  sinjT^-— z,'^z^z, 

Q'=0,  z>z,  or  z<Zc 

where  the  function  q{y),  which  specified  the  meridio- 
nal distribution  of  heating,  is  assumed  to  take  the 
form  of  (6).  Solution  (10)  will  be  evaluated  by  speci- 
fying the  vertical  scale  of  the  heating  Az=z,—Zc 
= 14  km,  Z)=  2.3X10“*  s'*  (damping  time  = 5 days) 
and  S,  = 3S,,  where  S,  and  St  are  the  static  stability 
for  the  stratosphere  and  troposphere,  respectively. 

Since  Mode  I is  less  trapped  in  the  vertical  we  e.x- 
pect  that  it  may  be  more  prominent  away  from  the 
source  region.  Using  (10a)  we  display  the  vertical 
wave  energy  flu.x  due  to  pressure  work,  q>'w',  at  Zt  as 
a function  of  vertical  wavelength  L,  and  zonal  wave- 
number  s(  = ka,  where  a is  the  radius  of  the  earth), 
in  Fig.  4.  The  cross-hatched  area  in  the  lower  part 
of  the  diagram  indicates  the  short  vertical  wavelength 
cutoff  required  by  (15).  The  maximum  response  occurs 
at  zonal  wavenumber  1-2  and  vertical  wavelength 
30-35  km.  As  expected,  this  result  is  similar  to  that 


>7 
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Frc.  3.  Imaginary  vertical  wavenumber  Xv  as  a function  of  vertical  wavelength  and  damping  time. 


obtained  by  Chans;  (1976)  except  that  the  vertical  Fig.  5 shows  this  quantity  as  a function  of  Lz  and  5. 
wavelength  is  slightly  longer  and  wavenumber  2 be-  Here  the  zonal  scale  selection  clearly  favors  the  lowest 
comes  comparable  in  importance  to  wavenumber  1.  zonal  wavenumber,  but  all  vertical  wavelengths 
This  vertical  wavelength  range,  according  to  Fig.  2,  beyond  the  short-wave  cutoff  have  about  the  same 
implies  a Doppler-shifted  phase  speed  ~50  m s“‘,  in  response.  Thus  the  most  prominent  viscous  mode  e.x- 
reasonable  agreement  with  the  observed  stratospheric  cited  should  be  of  zonal  wavenumber  1 and  vertical 
Kelvin  waves.  wavelength  >17  km.  For  vertical  wavelength  ~17  km 

For  the  highly  trapped  Mode  II  the  vertically-  the  Doppler-shifted  phase  speed  is  ~ 10  m s~‘  accord- 
inteerated tropospheric  kinetic  energy  is  considered  ing  to  Fig.  2.  These  values  agree  quite  well  with  the 
an  appropriate  measurement  for  the  response  spectrum.  40-50  day  oscillation  observed  by  Madden  and  Julian 


1 2 3 4 5 6 


Fig.  4.  Vertical  wave  energ>'  dux  due  to  pressure  work  at  the  tropopause  level, 
as  a function  of  vertical  wavelength  L,  and  zonal  wavenumber  i for  Mode  I. 
Units  are  arbitrary. 
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K.E.  Re*ponM  for  MODE  I 


Fig.  5.  Tropospheric  kinetic  energ>-  response  for  Mode  II  as  a function  of  vertical 
wavelength  and  zonal  wavenumber.  Units  are  arbitrary. 


(1972).  However,  if  the  forcing  function  has  a white 
noise  distribution  in  frequency  the  maximum  Mode  II 
response  is  still  somewhat  smaller  than  the  maximum 
Mode  I response.  The  absence  of  detectable  10-15 
day  Kelvin  waves  in  the  troposphere  therefore  indi- 
cates that  the  thermal  forcing  is  of  a red  noise  nature. 
However,  this  requirement  of  red  noise  distribution  is 
not  .so  much  for  the  explanation  of  the  stratospheric 
Kelvin  waves,  as  suggested  by  Holton  (1973);  rather, 
it  is  for  the  explanation  of  the  dominance  of  low- 
frequency  oscillations  in  the  troposphere.  The  10-15 
day  Kelvin  waves  in  this  case  are  of  relatively  higher 
frequencies  and  can  become  prominent  only  above  the 
troposphere  where  the  lower  frequency  viscous  mode 
is  greatly  attenuated. 

The  above  results  are  summarized  in  Table  2 within 
the  “non-resonant  selection”  category.  The  stationaiy 
waves  related  to  the  Walker  circulation  or  the  winter 
monsoon  may,  on  the  other  hand,  be  viewed  as  due 
to  resonant  excitation  by  the  stationary  heat  sources, 
because  the  forcing  function  has  a pre-selected  zero 
ground  phase  speed.  By  assuming  a zonal  mean  wind 


of  — 5 m 3~‘  in  the  source  region  the  Doppler-shifted 
phase  speed  of  this  resonance  is  5 m s~‘.  Referring 
back  to  Fig.  2,  response  is  possible  only  in  the  Mode  II 
regime  with  a vertical  wavelength  of  ~27  km  for  the 
case  of  5-day  damping  time.  This  vertical  wavelength 
agrees  much  better  with  the  observed  monsoon  or 
Walker  circulation  than  Webster’s  (1972)  inviscid 
calculation  of  ~ 1 km.  In  addition.  Fig.  5 still  indicates 
that  zonal  wavenumber  1 should  be  clearly  preferred 
over  shorter  waves  if  the  heat  source  has  a white 
noise  distribution  in  the  zonal  .scales.  In  the  real 
atmosphere  the  differential  heating  due  to  land-sea 
differences,  as  in  the  case  of  monsoon,  or  due  to  sea- 
surface  temperature  variations,  as  in  the  case  of 
Walker  circulation,  is  usually  confined  to  lower  zonal 
wavenumbers;  so  it  is  not  surprising  to  observe  the 
large  zonal  scale  of  these  circulations.  On  the  other 
hand,  Egger  (1976)  reported  that  the  response  to  the 
imp>osed  heating  anomalies  in  a linear  model  almost 
always  favors  the  lowest  wavenumbers  regardless  of 
the  zonal  scale  of  the  anomaly  impmsed  on  the  model. 
His  experimental  result  is  entirely  understandable  in 


Table  2.  The  most  eflSciently  excited  Kelvin  waves  predicted  by  the  viscous  theory. 


Vertical 

wavelength 

(km) 

Zonal 

wavenumber 

Doppler-shifted 
phase  speed 
(m  s-‘) 

Ground 
phase  speed 
(m  s"‘) 

Mode  I 

31-35 

1-2 

50-58 

45-53 

Resonant  selection 

Mode  U 

Mode  I 

^17 

(not  excitable) 

1 

15 

>10 

Non-resonant  selection 

Mode  II 

27 

1 

5 

0 

Damping  time  » 5 days. 

Vertical  scale  of  heating  • 14  km. 
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view  oi  our  result.  This  low -wavenumber  preference 
for  the  thermally  forced  Kelvin  wave  response  is 
related  to  (11)  and  the  explanation  has  been  given 
by  Webster  (1973)  and  Chang  (1976), 


damping  term,  i.e., 


4.  Further  comparison  with  the  40-50  day 
oscillation 

As  mentioned  previously,  the  two  solutions  of  (12) 
are  separated  by  the  curve  c,=  ei.  Since  for  Mode  I 
the  inertial  time  scale  is  faster  than  the  damping  time 
scale,  it  follows  that  the  main  balance  in  Eq.  (1)  is 
between  the  first  two  terms,  i.e.. 

In  this  case  u and  <t>  are  in-phase  and  the  wave  is 
characteristic  of  the  internal  gravity  waves.  For 
Mode  II  the  damping  time  scale  is  faster  so  that  the 
main  balance  is  between  the  last  two  terms,  i.e., 

0=  — ik4>—  Du. 

Thus  u and  (p  have  a quarter-cycle  phase  difference 
which  is  characteristic  of  a highly  viscous  wave  motion. 
If  we  e.tamine  the  schematic  diagram  of  the  40-50  day 
oscillation  (Fig.  1)  depicted  by  Madden  and  Julian, 
it  can  be  seen  that,  for  most  of  the  categories,  u and  <t> 
have  a phase  relationship  between  in-phase  and 
quarter-cycle  difference.  The  important  effect  of 
damping  is  thus  clearly  indicated. 

If  the  vertical  scale  is  large,  the  main  balance  in 
the  thermodynamic  energy  equation  is  between  the 
adiabatic  cooling  and  the  diabalic  heating.  The  residue 
is  approximately  balanced  by  the  advective  term  for 
Mode  I,  i.e., 

Q 

— ikcT=^ wT 

Cp 

and  T,  w will  have  a quarter-cycle  phase  difference. 
For  Mode  II  this  residue  is  mainly  balanced  by  the 


In  this  case  T,  w will  be  in-phase  if  diabatic  heating 
exceeds  adiabatic  cooling  and  half-cycle  out-of-phase 
if  adiabatic  cooling  dominates  over  diabatic  heating. 
In  Madden  and  Julian’s  (1972)  observational  study 
they  found  that  the  surface  pressure  is  half<ycle 
out-of-phase  with  the  300  mb  temperature  and  in-phase 
with  the  100  mb  temperature.  Since  vertical  motion 
is  close  to  in-phase  with  surface  pressure  as  shown 
by  Fig.  1,  it  may  be  inferred  that  T and  w are  in-phase 
at  300  mb  and  half<ycle  out-of-phase  at  100  mb.  The 
former  level  is  in  the  upper  troposphere  where  diabatic 
heating  must  dominate  for  the  motion  to  be  self- 
sustaining,  while  the  latter  level  is  near  the  tropopause 
where  the  diabatic  heating  vanishes  but  the  vertical 
motion  may  overshoot  through  the  cloud  top.  The 
thermal  structure  of  the  40-50  day  oscillation  thus 
also  appears  to  be  consistent  with  the  viscous  theoiy. 

When  damping  is  included,  the  meridional  structure 
of  the  Kelvin  waves  as  expressed  by  (6)  contains  a tilt 
of  the  wave  axis  because  t is  complex,  i.e.. 


Cr 

2 c,^+c. 


2Cr^+C^  J 


/ liXr  \ / d X,  Ci  \ 


(18) 


Here  the  first  part  of  the  right-hand  side  gives  a 
Gaussian  distribution  in  the  meridional  direction  with 
a Gaussian  folding  with  ya  given  by 


Fio.  6.  .Mean  phase  angles,  coherence-squares  and  background  coherence-squares  for  approximately  the  36-50  day  period  range  of 
cross  spectra  between  all  stations  and  Canton,  The  plotting  model  is  given  in  lower  right-hand  comer.  Positive  phase  angles  at  a station 
means  the  Canton  series  leads  that  of  the  sUtion,  Stations  indicated  by  a star  have  coherence-squares  above  the  background  at  the 
V5G?  level.  Mean  coherence -squares  at  Shemya  (52'’43'N,  174*’6'E)  and  Campbell  I.  (52®33'S,  IbP^Q'E)  Qnot  shown]  are  0,08  and  0,02, 
respectively.  Both  are  below  their  average  background  coherence-squares,  (From  Madden  and  Julian,  1972.) 
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Fig.  7.  Comparison  of  the  meridional  phase  tilt  of  the  40-50  day  oscillation  at 
140°E,  170°N,  observed  by  Madden  and  Julian  (1972),  and  the  value  of  the  most 
efficient  Mode  II  response  predicted  by  the  viscous  theory. 


This  meridional  trapping  scale  is  identical  for  both 
the  regular  and  the  viscous  modes  for  a given  vertical 
wavelength.  In  the  case  of  the  most  efficient  Mode  II 
response  discussed  in  the  previous  section,  \yo'. 
= 1500  km  and  agrees  quite  well  with  the  observed 
40-50  day  oscillation  The  second  part  of  the  right- 
hand  side  of  (18)  gives  rise  to  a periodic  variation  of 
the  waves  in  y,  and  therefore  a meridional  tilt  of  the 
wave  axis.  It  also  implies  a poleward  phase  propaga- 
tion for  the  Kelvin  waves  because  Cr>0,  c,>0  and 
X,>0. 

A strong  meridional  tilt  was  found  by  Madden  and 
Julian  for  the  40-50  day  oscillation  using  cross-spectral 
analysis  of  the  surface  pressure  data.  Their  phase 
diagram,  constructed  by  using  the  Canton  Island  series 
as  the  base  series,  is  reproduced  here  in  Fig.  6.  Al- 
though the  meridional  phase  tilt  seems  to  vary  in 
different  geographical  locations,  near  Canton  Island 
and  near  the  equator  the  tilt  is  primarily  NVV-SE 
north  of  the  equator  and  NE-SW  south  of  the  equator. 
For  eastward  propagating  waves  this  indicates  a pole- 
ward  phase  propagation.  In  Fig.  7 the  phase  tilt  of 
the  40-50  day  oscillation  deduced  from  Fig.  6 at  two 
different  longitudes,  170°\V  (near  Canton  Island)  and 
140®E  (active  convection  area  in  Indonesia),  are 
plotted  along  with  that  given  by  the  viscous  theory. 
Within  the  limitation  of  our  simple  theory  and  the 
approximations  invoked,  the  theoretical  result  can  be 
considered  as  fairly  good  in  simulating  the  observed 
phase  tilts. 


5.  Concluding  remarks 

We  have  shown  that  the  inclusion  of  a damping 
term  can  result  in  two  types  of  dispersive  relationships 
for  the  forced  Kelvin  waves.  The  first  typ)e  (Mode  I) 
is  characteristic  of  the  regular  internal  gravity  waves 
which  propagate  rapidly  in  the  zonal  direction  and 
attenuate  slowly  in  the  vertical.  The  second  type 
(Mode  II)  is  dominated  by  the  viscous  damping  time 
scale  and  has  slower  zonal  phase  speed  and  faster 
vertical  attenuation.  This  viscous  mode  reduces  to  a 
trivial  solution  in  the  inviscid  limit  so  its  existence 
is  supported  only  by  damping.  Assuming  that  the 
main  thermal  sources  in  the  tropics  are  confined  to 
the  troposphere  and  have  a red  noise  distribution 
which  enhances  the  response  of  the  viscous  mode, 
this  theory  predicts  that  the  stratospheric  waves 
behave  like  the  regular  mode  and  the  tropospheric 
waves  behave  like  the  viscous  mode.  Our  theory  in- 
cludes features  that  are  consistent  with  several  results 
found  by  previous  theoretical  or  numerical  studies  of 
the  equatorial  waves  that  incorporate  the  damping 
mechanism.  Lindzen  (1971)  and  Holton  and  Lindzen 
(1972)  found  that  as  a wave  is  propagating  in  weak 
vertical  wind  shear,  the  efficiency  of  damping  would 
change  inversely  according  to  the  change  of  the 
Doppler-shifted  phase  speed.  This  is  consistent  with 
our  result  that  the  slower  moving  waves  are  more 
trapped  in  the  vertical.  In  a numerical  diagnostic 
model  of  the  Kelvin  waves,  Holton  (1973)  observed 
a maximum  response  in  the  troprospheric  wavenum- 
ber 1 sprectrum  with  preriod  >30  days.  This  response 
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seems  to  correspond  to  our  viscous  mode.  Very  recently 
Stevens  el  d.  (1977)  reported  that  with  a parameter- 
ization of  the  cumulus  damping,  the  tropical  waves 
have  very  little  vertical  structure.  This  is  again  ex- 
plainable by  our  theory  because  the  vertical  wave- 
length is  increased  as  a result  of  damping. 

In  real  atmosphere  the  damping  due  to  cumulus 
convection  and  other  processes  is  much  more  com- 
plicated than  the  simple  linear  drag  term  used  in  our 
analysis.  However,  the  agreement  with  observed  zonal 
wind  oscillations  such  as  the  stratospheric  Kelvin 
waves,  the  40-50  day  oscillations,  and  the  stationary 
monsoon  and  Walker  circulations  strongly  indicate 
that  our  theoretical  interpretation  is  quite  relevant 
to  these  oscillations.  In  addition,  although  the  present 
analysis  is  limited  to  Kelvin  waves  only,  the  vertical 
structure  equation  applies  to  all  equatorial  waves. 
Hence  it  would  be  quite  possible  that  the  general 
result  also  applies  to  other  thermally  controlled  tropical 
wave  disturbances.  This  is  particularly  interesting  in 
view  of  the  large  vertical  scale  and  slow  Doppler- 
shifted  phase  speed  usually  observ'ed  in  the  tropical 
tropospheric  waves  such  as  the  easterly  waves.  In  a 
discussion  of  the  present  results,  Wallace  (personal 
communication)  suggested  that  they  would  also  ex- 
plain the  seemingly  non-dispersive  westward  propaga- 
tion of  the  synoptic-scale  cloud  clusters  observed  in 
the  time  sections  of  tropical  satellite  pictures  (Chang, 
1970).  In  these  sections  the  propagating  cloud  patterns 
have  quite  uniform  phase  speeds  but  the  apparent 
wavelength  varies  over  a large  range.  This  non-dis- 
persive  character  led  Wallace  (1972)  to  propose  that, 
despite  the  strong  vertical  and  longitudinal  variations 
of  the  basic  flow,  they  are  more  likely  a result  of 
advection  by  basic  flow  than  a manifestation  of  some 
kind  of  Rossby  waves.  He  argued  that  the  difficulty 
of  ascertaining  a steering  level  is  not  as  severe  as  that 
of  explaining  why  the  beta-effect  does  not  produce 
dispersive  phase  speeds.  However,  to  interpret  these 
patterns  as  being  passive  cloud  masses  advected  by 
the  basic  flow  is  also  unsatisfactory  due  to  the  fact 
that  they  are  frequently  observed  to  be  associated 
with  active  convections  and  that  they  sometimes 
decay  and  regenerate  during  the  course  of  propaga- 
tion. These  difficulties  are  partially  alleviated  in  the 
framework  of  the  present  tlisory  tecause  the  viscous 
mode  should  have  a much  smaller  Doppler-shifted 
phase  speed  and  would  appear  as  if  they  are  advected 
by  the  basic  easterly  current  at  some  levels. 
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Some  Theoretical  Problems  of  the 
Planetary-Scale  Monsoons 

By  C iim-Pia  C manc;'  ) 

Sumtnar\  St>mc  iinporiani  theoretical  problems  ol  the  planeiarv-scale  monsoons  v^hich  have  arisen 
from  recent  advances  ol  observational  studies  are  reviewed  These  include.  ( 1 ) the  requirement  ol'a  strong 
damping  mechanism  m the  planetary  scale  vorlicily  budeel  oCsumnier  monsoon  and  a similar  hut  weaker 
requirement  lor  the  winter  nn'iisoon.  (2)  the  localized  barotropic  instabiliiv  of  the  summer  monsi)on 
which  IS  a result  ol  the  strong  zonal  asymmetry  ofthe  planetary-scale  flow  and  causes  signiiicani  nonlinear 
enerev  conversnms:  and  ( 'o  the  oscillations  o|‘  the  planetary-scale  monsoons.  It  is  pointed  out  that  these 
problems  are  inier-relaied  and  their  understandine  is  also  important  tor  the  proper  simulation  of  other 
scales  of  motion  ol  the  monsinm  circulation. 

key  words;  Monsoon  (.)scillati(m.  \orticiiy  budget;  Barotropic  instabiluy 


/.  liUriHha  titiii 

The  monsoon  circulation  of  both  the  northern  summer  and  the  northern  vvmter 
(hereafter  the  word  'northern'  will  be  neglected)  are  characterized  by  large  planetary 
scale  quasi-stationary  waves  in  the  subtropics  and  tropics.  .As  have  been  shown  by 
KKisuNAMi  Rri  (I9"lai.  KRisuNvMt  rti.  Kana.mitsi  . Koss  and  Li;i;  (I9'^?b)  and 
others,  these  waves  are  a result  of  the  thermally  direct  overturnings  induced  by  the 
dilferential  heating  associated  with  the  land-sea  distribution.  The  major  component 
of  the  heating  is  the  convective  latent  heat  rclea.se  which  is  most  intense  over  the 
Tibetan  and  Mexican  highlands  during  summer  and  the  maritime  continent  of 
the  Malaysia-lndonesia  region  during  winter.  .Mainly  due  to  the  availability  of  the 
commercial  aircraft  data,  the  principal  features  of  these  zonal  asymmetries  have  been 
best  depicted  by  the  w ind  analyses  near  2(X)  mb  by  KRisuNAMURTiand  Roghrsi  1970). 
KRisHNAMLRri  c/  uL  ( 1973b)  and  Sadi.kr  (1975).  During  summer  these  features 
include  a strong  anticyclone  centered  over  the  Tibetan  plateau,  a weaker  anticyclone 
over  Mexico  and  elongated  troughs  oriented  northeast  southwest  in  the  mid-Pacific 
and  mid-Atlantic.  These  are  all  well  presented  in  Fig.  1 which  shows  the  time-mean 
streamfunction  at  2(X)  mb  for  June  August  1967  as  determined  by  Krishnaml  rti 
( I9tla).  During  winter  a large  streamfunction  high  over  the  equatorial  South  China 
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Figure  I 

Summer  mean  (June  August  PM7)  2(MFmb  iireamiunction,  contour  interval  is  50  ^ 10*  m*  sec  ‘.(From 

KRISHNA.MLRn.  I97la.) 


to » 0 50 40 90 120  ISO 


go  ao 


«0  50 


50«0  90  f20iS0«0l90a0 


I’lgure  2 

Winier  mean  ( December  February  I9?h  20tumb  sireamlunction.  contour  interval  is  50  10*  m- sec  ' 

(From  Krishnamlrm  ef  u/.  19‘’3b.) 

Sea  and  western  Pacilie  replaces  the  Tibetan  high  to  become  the  most  pronounced 
planetary  scale  feature,  with  two  weaker  highs  appearing  over  the  equatorial  eastern 
Africa  and  the  northeastern  South  America  and  three  oceanic  troughs  extending 
from  the  equator  ro  the  southern  subtropics  These  features  are  shown  in  Fig.  2 which 
is  the  200  mb  time-mean  streamfunction  for  the  winter  of  1971  determined  by 
Krishnamlrti  <■/  ill.  ( 1973b). 

A number  of  interesting  theoretical  problems  have  arisen  from  the  observed 
mon.soon  circulations.  In  this  paper  we  will  review  some  of  the  problems  that  are 
related  to  the  planetary-scale  monsoon.  They  include; 


( 1 1 The  requirement  of  a strong  damping  mechanism  in  the  planetary  scale 
vorticity  budget. 
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(2)  The  localized  barotropic  instability  associated  with  the  pliinclary-sctilc 
summer  monsoon,  and 

(3)  Oscillations  of  the  planetary-scale  monsoons. 


2.  The  reiiuirement  oj  a strong  planetary  scale  clamping  mechanism 

In  a linear,  one-level  numerical  model.  Hoi.ton  and  Coi.jon  (1972)  used  the 
barotropic  vorticity  equation  to  diagnose  the  budget  of  the  seasonal-mean  planetary 
scale  vorticity  at  200  mb  during  summer.  They  used  the  200  mb  mean  divergence 
observed  by  Krishnami  rti  (1971a.  b)  as  the  steady  forcing  function  with  a mean 
zonal  flow  also  specified  by  the  observed  values.  Figure  3 shows  the  velocity  potential 


[■Igure  5 

Summer  mean  (June  August  1%7)  2ll0-mb  velocity  potential,  contour  interval  is  It)  « It)'  m-  sec  ' The 
lew  streamlines  sketched  indicate  the  direction  of  the  mean  divergence  motions.  (From  Krishnami  nri. 

|97|a.i 


computed  by  Krishnamurti  (1971a)  which  expresses  the  observed  mean  divergence 
field.  Comparing  this  figure  with  Fig.  1.  it  is  evident  that  the  time-mean  vorticity  and 
divergence  fields  are  approximately  half-cycle  out  of  phase.  Flowever.  the  resultant 
steady-state  planetary  scale  vorticity  field  calculated  by  Holton  and  Colton  is 
about  one-quarter  cycle  out  of  phase  from  the  divergence  field,  unless  a very  strong 
damping  term  is  added  to  the  mean  vorticity  equation.  This  difficulty  arises  because, 
in  the  absence  of  strong  damping,  the  main  balance  in  the  vorticity  equation  is  be- 
tween the  generation  by  horizontal  divergence  and  the  advection  by  mean  flow.  This 
near-balance  requires  that  the  vorticity  be  one-quarter  cycle  out  of  phase  from  the 
divergence.  The  very  efficient  generation  by  the  strong  divergence  also  causes  the 
amplitudes  of  the  long  waves  to  be  significantly  larger  than  those  observed.  Figure  4 
(from  Holton  and  Colton,  1972)  compares  the  diagnosed  and  observed  amplitudes 
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and  phases  of  the  streamfunction  for  zonal  wavenumber  one  as  a function  of  latitude. 
When  the  linear  damping  coefficient  is  increased  from  a modestly  strong  value  of 
5 X lO'*’  sec*  ‘ to  a very  strong  value  of  1.5  x 10'*  sec“‘,  the  discrepancy  dis- 
appears and  the  phase  and  amplitude  of  the  long  waves  become  in  close  agreement 
with  observations. 


PHASE  (LONGITUOe) 

9or«  O'  ao'E  leo"  sew 


amplitude 

Figure  4 

Amplitude  in  units  of  Itr  m’  sec  ‘ (left-hand  curves)  and  phase  (right-hand  curves)  of  the  streamfunction 
for  zonal  wavenumber  I as  a function  of  latitude.  Solid  curve  is  observed  streamfunction;  short  and  long 
dashes  are  the  solutions  computed  by  Holton  and  Colton  ( 19^2)  for  damping  rate  of  1.5  ' 10  ’sec  ' 
and  5 ' lO  " sec  ',  respectively  (From  Holton  and  Colton.  1972.) 


Hijlton  and  Colton's  (1972)  findings  suggest  that  a very  strong  sink  with  a 
damping  time  scale  faster  than  1 day  is  required  to  balance  the  linear  vorticity  budget 
of  the  seasonal-mean  200  mb  planetary-scale  flow  auring  summer.  This  damping  is 
much  stronger  than  the  normal  dissipation  process  in  the  upper  troposphere  whose 
damping  time  is  on  the  order  of  weeks.  In  looking  for  an  explanation  of  this  excessive 
damping  requirement,  they  suggested  that  two  possible  mechanisms  not  included  in 
their  model  may  be  responsible  for  the  strong  damping.  The  first  is  the  vertical  mixing 
of  vorticity  by  cumulonimbus  convection  which  is  very  strong  over  the  Asian  conti- 
nent and  Mexican  highland  during  the  northern  summer.  This  mechanism  has  been 
suggested  by  Riehl  and  Pearce  (1968).  Reed  and  Recker  (1971).  Wallace  ( 1971). 
WiLLiA.vis  and  Gray  (1973)  and  others  to  balance  the  vorticity  budget  of  the  synoptic- 
scale  easterly  waves  in  the  tropical  troposphere.  However.  Holton  and  Colton 
(1972)  felt  that  it  alone  is  still  insufficient  to  account  for  the  required  fast  damping. 
They  then  proposed  the  second  mechanism  which  is  related  to  the  day-to-day 
fluctuations  of  the  monsoon  vorticity  and  divergence  fields  observed  by  Krishna.mcr- 
Ti  (1971b).  Considering  deviations  from  the  time  mean,  the  complete  divergence  term 
in  the  time-mean  vorticity  equation  is  <(l  -f-  y)(5>  = (<C>  + y)<(5>  ^ <'(5  C’>  where  ^ 
and  / are  the  relative  and  planetary  vorticities.  respectively;  <5  is  divergence  and  the 
primes  are  deviation  from  the  time  mean  < >.  The  last  term,  which  was  excluded  in 
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Holton  and  Colton's  linear  calculations,  tends  to  be  negative  according  to  the 
observed  out-of-phase  correlation  of  the  fluctuating  vorticity  and  divergence  fields. 
This  term  therefore  may  counteract  the  vorticity  generation  by  the  time-mean 
divergence.  Holton  and  Colton  then  speculated  that  the  combination  of  this 
fluctuation  correlation  and  the  cumulus  mixing  mechanisms  could  be  responsible  for 
the  less  than  one  day  damping  time. 

The  requirement  of  a strong  vorticity  sink  has  also  been  confirmed  by  other 
numerical  investigations  that  do  not  include  the  strong  damping.  As  a follow  up  of 
the  work  of  Holton  and  Colton  (1972).  Colton  (1973)  used  a semi-spectral  repre- 
sentation of  the  one-level  nonlinear  barotropic  vorticity  equation  to  study  the  scale 
interactions  of  the  200  mb  summer  monsoon  flow  forced  by  the  observed  mean 
divergence.  .As  a consequence  of  the  inclusion  of  only  modest  damping  (time  scale 
5 days),  the  resultant  planetary-scale  waves  have  amplitudes  considerably  larger  than 
the  seasonal  mean  values.  Westward  phase  shifts  of  50°  for  the  Tibetan  high  and  30° 
for  the  mid-Pacific  trough,  in  comparison  with  the  observed  positions,  are  aLso 
produced  by  the  model  (Fig.  5).  .Abboit  ( 1973)  developed  a 3-level  nonlinear  quasi- 
geostrophic  hemispheric  model  to  simulate  the  time-mean  climatology  of  the  summer 


Figure  5 

Streamfunction  of  the  forced  motion  averaged  over  days  15  through  40  of  the  numerical  integration  by 
flit. niN  ( 1977).  Contour  interval  IS  50  ' lO' m’ sec 


monsoon.  Todrive  the  monsoon  circulation,  he  specified  a steady  heating  distribution 
which  was  found  by  inverting  the  governing  equations  with  the  vorticity  field  giving 
by  the  observed  time-mean  values.  This  heating  distribution  (Fig.  6)  in  turn  success- 
fully simulates  the  principal  features  of  the  planetary-scale  flow  as  well  as  the  nonlinear 
energy  exchanges.  However,  it  is  about  a quarter  cycle  out  of  pha.se  from  that  implied 
by  the  observed  mean  divergence  field  and  the  satellite  cloud  data.  This  apparently  is 
also  due  to  the  absence  of  a strong  damping  mechanism  in  the  model  w hich  forces  a 
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Figure  6 

Time-independent  horizontal  distribution  of  the  diabatic  heating  rate  at  35(i  mb  as  computed  by  AaBor 
( 197.1).  Isopleths  arc  at  5 < 10  ^ m'*  sec  '. 

balance  between  the  advection  and  generation  of  vorticity.  In  a 3-component  semi- 
spectral  primitive  equation  model.  Muraka.vii  (1974)  also  used  a specified  heating 
function  to  study  the  planetary  scale  monsoon  flow  during  summer.  Without  the 
inclusion  of  the  strong  damping  and  allowing  the  specified  heating  to  fluctuate  with 
time,  he  was  able  to  obtain  the  proper  amplitude  and  phase  of  the  long  waves. 
However,  the  resultant  divergence  field  in  his  model  is  somewhat  different  from  that 
observed  by  Krishnavu  rti  (1971a).  The  greatest  discrepancy  occurs  at  200  mb. 
where  Krishnamurti  (1971a)  found  very  strong  divergence  in  the  Tibetan  high  area 
during  the  1967  summer,  while  in  Murakami's  ( 1974)  model  this  level  is  almost  non- 
divergent.  The  model  also  produces  a different  horizontal  distribution  of  long  wave 
divergence  compared  to  the  observed  1967  distribution,  with  the  Tibetan  divergence 
center  sh'fted  slightly  to  the  east  and  the  mid-Atlantic  convergence  center  shifted  to 
the  longitude  of  central  Africa.  Thus  the  magnitude  of  divergence  forcing  is  much 
reduced  and  a phase  shift  between  divergence  and  anticyclonic  vorticity  also  results  ’ ). 
In  addition,  the  advective  effect  is  not  properly  treated  because  the  zonal  mean  flow 
is  not  allowed  to  interact  with  the  three  interacting  long  waves.  Hence  it  is  unclear 
what  is  the  main  mechanism  that  produces  the  correct  amplitude  and  phase  of  the 
vorticity  field. 

Looking  for  alternative  sources  for  the  strong  damping,  Barcilon  and  Ccxjper 
(1974)  hypothesized  that  a thermal  boundary  layer  effect  could  operate  in  the 
atmosphere  during  summer  monsoon.  In  their  simple  analytical  model  a sinu.soidal 
temperature  variation  is  applied  at  the  lower  boundai’y  which  provides  the  planetary 
scale  differential  heating.  A number  of  solutions  were  found  for  a balance  between 


•I  However.  Mi  rakami  (personal  communication)  found  thal  this  honzonlal  distribution  ofdiser- 
zence  is  in  good  agreement  with  the  1970  summer  data. 
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the  horizontal  advection  and  vertical  diffusion  of  heat  in  the  thermodynamic  energy 
equation,  and  one  of  them  allows  a half-cycle  out  of  phase  shift  between  the  upper 
level  vorticity  and  divergence.  However,  since  the  monsoonal  differential  heating  in 
the  real  atmosphere  is  mainly  due  to  the  latent  heat  release  of  cumulus  convection 
which  is  maximum  in  the  middle-upper  troposphere,  it  is  uncertain  whether  Barcilos 
and  Cooper's  (1974)  theory  is  completely  applicable  to  the  real  atmosphere. 

Another  explanation  for  the  observed  vorticity -divergence  phase  relationship  was 
offered  by  Paegle  and  Paegle  (1976).  They  speculated  that  a high  pressure  region 
with  strong  pressure  gradients  may  cause  the  geostrophic  vorticity  to  become  less 
than  one-half  of  the  negative  planetary  vorticity.  This  implies  that  the  balance 
equation  becomes  non-ellipticand.  under  certain  rather  strong  assumptions,  leads  to  a 
rapid  vorticity  adjustment  with  strong  horizontal  divergence  occurring  in  the  negative 
relative  vorticity  region.  However,  the  resultant  steady-state  absolute  vorticity 
vanishes  and  this  ethciently  eliminates  the  generation  of  vorticity  by  divergence.  Such 
large  negative  v alues  of  relative  v orticity  are  not  observed  in  the  large-scale  monsoon 
flow.  In  fact  the  minimum  seasonal  mean  relative  vorticity  at  200  mb  during  the  1967 
summer  as  observed  by  Krisunavu  rti  is  ^ ■<  10"  '' sec  ‘ which  occurs  in  the 

vicinity  of  the  Tibetan  and  Mexican  highs  near  .10°N.  This  magnitude  is  considerably 
smaller  than  that  of  the  Coriolis  parameter.  Thus  it  appears  unlikely  that  this  super- 
gradient effect  is  relevant  to  the  planetary-scale  monsoon  problem. 

To  test  the  hypotheses  proposed  by  Holton  and  C(jlton  (1972).  Chang  and 
Pi.NTiMONTi  (1977)  used  a .'^-level  primitive  equation  model  with  a 4°  tinite-dilference 
resolution  in  both  the  zonal  and  meridional  directions.  The  zonal  mean  (low  is 
included  usinga  Newtonian  term  which  maintains  the  proper  north-south  temperature 
gradient  by  continuously  adjusting  each  temperature  toward  a thermal  equilibrium 
determined  by  the  zonaily  symmetric  climatology.  The  zonally  asymmetric  motions 
are  generated  by  a specified  perturbation  heating  function  whose  horizontal  distribu- 
tion follows  that  of  the  200  mb  seasonal  mean  divergence  observed  by  Krishnaml  rit 
(19"lb).  It  turns  out  that  this  heating  distribution  almost  reproduces  the  observed 
mean  divergence  in  the  upper  troposphere,  which  suggests  that  the  choice  of  the 
forcing  function  is  proper.  Two  experiments  were  carried  out  by  integrating  the  model 
in  time  until  25  days  after  it  reached  the  respective  steady  states.  Both  experiments 
include  only  weak  damping  (time  scale  >7  days)  in  the  free  atmosphere  but  the 
heating  function  was  held  steady  only  in  the  first  experiment  and  allowed  to  fluctuate 
in  time  with  an  amplitude  equal  to  one-half  of  the  time-mean  heating  in  the  second 
experiment.  The  fluctuation  period  was  set  to  10  days.  Comparing  the  results  of  the 
two  experiments.  Chang  and  Pentimonti  found  essentially  no  difference  between 
the  time-mean  planetary  scale  fields  except  a slightly  reduced  amplitude  for  the 
fluctuating  case.  Figure  7 is  their  time-mean  w ind  field  at  250  mb  (near  the  upper  level 
of  the  model)  for  the  fluctuating  heating  case.  It  can  be  seen  that  the  simulated 
Tibetan  and  Mexican  anticyclones  are  shifted  about  .10°  and  25°.  respectively,  from 
the  observed  positions.  These  phase  shifts  are  almost  identical  to  those  produced  by 
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their  steady-heating  case.  Thus  the  fluctuating  component  as  specified  in  the  model 
is  almost  inconsequential  in  providing  the  damping  needed  for  the  phase  shift  prob- 
lem. On  the  other  hand,  in  both  e.xperiments  the  amplitudes  of  the  long  waves  are  not 
e.xcessively  large  compared  to  observations.  This  result  is  similar  to  that  of  Abbott 
(1973)  and  seems  to  suggest  that  nonlinear  proces.ses  could  be  responsible  for  part  of 
the  planetary-scale  vorticity  sink  deduced  by  Hoi.TosandCot.TON'slinearcalculations. 
The  fluctuation  of  the  planetary-scale  flow  apparently  enhances  these  nonlinear 
conversions  resulting  in  the  further  reduction  of  its  time-mean  amplitude. 

A noteworthy  feature  in  C'hanci  and  PiiM  i.sK inti's  ( 1977)  time-mean  fields  is  that 
the  positions  of  the  upper  tropospheric  mid-oceanic  troughs  do  not  have  significant 
phase  errors  as  those  of  the  anticyclones.  As  pointed  out  by  Kkisiinami  rti.  DACitii  - 
PATV.  Fits.  KASAMnst  and  Lfct  (1973a).  the  damping  due  to  cumulus  transports 
should  not  operate  over  the  relatively  cloud-free  oceanic  areas.  With  the  inclusion  of 
nonlinear  processes  (and  also  large-scale  vertical  advections).  the  damping  require- 
ment implied  by  the  phase  errors  of  Chanc;  and  Piati.monti's  time-mean  upper  level 
anticyclones  is  weaker  than  that  found  by  Hotton  and  Colton  (I9'?2).  Perhaps  a 
damping  time  longer  than  1 day.  which  seems  more  reasonable  than  Holton  and 
f dlion's  value  if  interpreted  as  due  to  the  cumulus  transports,  may  be  suflicient. 
In  this  case  a numerical  experiment  similar  to  C hano  and  Pi  ntimonli's.  hut  in- 
corporating a moderately  strong  damping  term  with  time  scale  order  of  a few  days  in 
the  continental  convective  area  only,  may  provide  some  insight  as  to  whether  cumulus 
transport  processes  are  responsible  tor  the  missing  vorticity  sink. 

The  studies  mentioned  above  all  concern  the  summer  monsoon  only.  During 
winter  Krishnami  rti  ct  <;/.  (I9~3b)  observed  that  the  rotational  component  of  the 
seasonal-mean  flow  of  the  active  monsoon  region  is  almost  one  order  smaller  than 
that  of  the  summer  monsoon,  while  the  divergent  component  has  a comparable 
magnitude.  The  reason  for  this  appears  to  be  the  proximity  of  the  winter  monsoon 
convections  to  the  equator,  which  reduces  the  etliciencv  of  vorticity  generation  by 
horizontal  divergence.  However,  the  importance  of  a dissipation  mechanism  in  the 
winter  monsoon  vvas  also  brought  out  m a few  theoretical  studies  recently.  In  a linear 
study  of  the  stationary  forcing  in  the  tropical  atmosphere.  Wlbstlr  (I9"’2.  |973) 
found  that  the  primary  response  in  the  equatorial  region  corresponds  to  the  atmos- 
pheric Kelvin  waves  (Hoi.roN  and  Lind/.i  n.  196S).  .A  theoretical  calculation  shows 
that  the  zonal  wavelength  of  this  Kelvin  wave  respon.se  is  very  long  while  its  vertical 
wavelength  is  very  short  [0(1  km)].  The  former  agrees  quite  well  with  observed 
equatorial  stationary  circulations  such  as  the  winter  monsoon,  but  the  latter  is  one 
order  of  magnitude  shorter  than  the  observed  vertical  wavelengths.  This  discrepancy 
arises  because  the  equatorial  Kelvin  waves  behave  like  internal  gravity  waves  in  the 
zonal-vertical  plane  and  their  phase  speed.  < . is  related  to  an  equivalent  depth  scale  h 
by  the  relationship  < = (<//i)'  ’.  These  stationary  waves  are  imbedded  in  a mean 
easterly  zonal  current  so  that  they  have  a Doppler-shifted  phase  speed  of  a few  meters 
per  second,  w hich  implies  a very  short  vertical  wavelength.  To  get  around  this  problem 
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in  simulating  the  observed  stationary  motions.  Webster  (1972,  1973)  adopted  a two 
level  numerical  model  with  friction  which  is  quite  successful  in  reproducing  many  of 
the  observed  features. 

Very  recently  Chang  (1977)  re-examined  this  problem  by  incorporating  linear 
damping  terms  in  the  thermally  forced  Kelvin  wave  equations.  The  dispersive  rela- 
tionship of  the  Kelvin  waves  may  be  written  as  </  = 5‘  where  /.  is  the  vertical 
wavenumber  and  5 is  a static  stability  parameter.  The  linear  damping  introduces  an 
imaginary  part  to  c which  causes  /.  to  become  complex  also.  Consequently  the 
dispersive  relationship  becomes  quadratic  and  has  two  types  of  solutions.  The  first 
type  arises  when  the  inertial  time  scale  is  faster  than  the  damping  time  scale.  This 
solution  resembles  the  regular  internal  gravity  waves  whose  vertical  wavelength 
decreases  as  the  Doppler-shifted  phase  speed  decreases.  The  second  type,  called  the 
viscous  mode  by  Chang  (1977),  arises  when  the  damping  time  scale  dominates  over 
the  inertial  time  scale  and  its  vertical  wavelength  increases  as  the  Doppler  phase  speed 
decreases  (Fig.  8).  This  solution  is  also  much  more  trapped  in  the  vertical  so  that  it  has 
appreciable  amplitude  only  in  the  region  of  thermal  forcing  (troposphere).  In  the 
in  viscid  limit  this  solution  reduces  to  a trivial  solution  with  vanishing  amplitude.  Fora 
reasonable  range  of  physical  parameters.  Chang  ( 1977)  calculated  a vertical  wave- 
length of  >17  kms  for  the  stationary  response,  which  is  a significant  improvement 
over  the  inviscid  theory.  The  horizontal  structure  of  this  solution  is  also  in  good 


Figure  S 

Real  Doppler-shified  phase  speed.  ( as  a runction  of  vertical  wavelength.  and  damping  time  lor  the 

Kelvin  waves.  (From  Chanu,  1977.) 
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agreement  with  observations  and  Webster's  (1972,  1973)  numerical  model.  This 
requirement  of  a damping  mechanism  with  time  scale  on  the  order  of  5 days  is  much 
weaker  than  that  found  by  Holton  and  Colton  (1972)  for  the  summer  monsoon 
and  could  well  be  e,xplained  by  cumulus  transports.  Detailed  observational  studies  of 
the  winter  monsoon  are  needed  to  improve  the  theoretical  understanding  of  this 
problem. 


3.  Localized  harotropic  instahiUty  of  the  planetary  scale  summer  monsoon 

As  mentioned  previously,  the  amplitude  of  stationary  long  waves  during  summer 
is  better  simulated  by  the  nonlinear  models  than  by  the  linear  models  when  excessively 
large  damping  is  excluded.  This  apparently  is  due  to  the  barotropic  energy  conv  ersions 
from  planetary  to  both  zonal  and  synoptic  scales  which  are  allowed  in  the  nonlinear 
models  only.  Observational  evidence  of  this  type  of  energy  conversion  has  been 
presented  by  Kanamitsu.  Krishnaml  rti  and  Depradine  (1972).  who  computed  the 
wave-zonal  flow  and  wave-wave  energy  exchanges  for  the  tropical  belt  from  1 5°S  to 
I5°N  using  the  1967  summer  2(K)  mb  wind  data  analyzed  by  Krishnaml  rti  and 
Rogers  1 1970).  In  particular,  they  showed  that  the  synoptic  scale  waves  receive  kinetic 
energy  at  a fairly  significant  rate  (order  of  10 m-  sec  •')  from  both  the  zonal  flow 
and  the  planetary-scale  waves.  .A  similar  result  for  the  equator  to  30“N  belt  was  also 
obtained  by  Krlshnaml  rti  ei  al.  ( 1973b). 

Transient  synoptic-scale  wave  disturbances  have  been  observed  in  the  tropical 
upper  troposphere  in  all  seasons.  During  summer,  the  most  active  season,  the  distur- 
bances are  found  to  occur  primarily  in  two  types  of  regions.  One  type  is  in  the  vicinity 
of  the  mid-oceanic  troughs  (Rick.  1959;  Sadler.  1967;  Miller  and  Carlson.  1970; 
Frank.  1970)  where  the  disturbances  may  assume  the  form  of  open  waves  or  closed 
vortices.  Occasionally  these  disturbances  may  develop  in  the  vertical  or  be  coupled 
with  low  level  waves,  and  intensification  leading  to  tropical  cyclone  development 
sometimes  may  also  follow.  .Another  region  of  transient  wave  activity  is  along  the 
upper  tropospheric  easterly  jet  over  the  Indian  ocean,  where  Krishna.mlrti  and 
Rcx3ers(  1970)  have  reported  fairly  regular  occurrences  of  synoptic-scale  disturbances. 
.Apparently  these  regions  are  where  most  of  the  barotropic  energy  conversion  from 
zonal  and  planetary  scales  to  the  synoptic  scale  takes  place.  Thus  it  appears  that  the 
nonlinear  energy  transfers  of  the  summer  monsoon  have  a persistently  localized 
nature.  Another  evidence  of  this  is  related  to  the  linear  barotropic  instability  criterion, 
which  has  been  shown  by  MERtLEES  (1968)  to  be  usually  satisfied  if  nonlinear  energy 
exchanges  are  very  prominent.  The  upper  tropospheric  zonal  mean  flow  of  the 
summer  monsoon  does  not  satisfy  this  criterion,  but  regions  of  large  negative 
gradients  of  relative  vorticity  (Fig.  9)  where  this  instability  criterion  is  satisfied 
locally  can  be  found  in  the  seasonal  mean  2(X)  mb  analysis  (Krishnamurti.  1971b). 
These  negative  vorticity  gradients  occur  mainly  in  the  vicinity  of  the  easterly  jet  and 


Figure  9 

Summer  mean  (June-August  1967)  relative  vorticity  at  2(X)  mb.  contour  interval  is  10  ’ sec  (From 

KRISHNA.MLRri,  1971b.) 


the  mid-oceanic  troughs,  and  on  some  individual  days  their  magnitude  may  even 
greatly  exceed  the  seasonal  mean  values.  Therefore  it  is  plausible  that  at  least  some  of 
the  transient  synoptic  wave  activities  in  these  regions  is  generated  by  the  localized 
(and  sometimes  short  term)  barotropic  instability,  which  is  clearly  a result  of  the 
zonal  asymmetry  of  the  planetary-scale  summer  monsoon. 

The  kinetic  energy  conversions  from  the  planetary  scale  and  the  zonal  flow  to  the 
synoptic  scale  in  the  upper  troposphere  have  been  simulated  numerically  by  Abbott 
(1973)  and  Krishnamlrti  et  al.  (1973a).  The  conversions  produced  by  Abbott's 
(1973)  model  are  roughly  of  the  same  order  of  magnitude  as  thos?  computed  by 
Kanamitsu  et  al.  ( 1972).  However,  the  conversions  obtained  by  Krishna.murti  ct  al. 
(1973a)  are  about  one  order  smaller.  This  appears  to  be  related  to  their  treatment  of 
the  Tibetan  high  as  a (ixed  barrier  which  prevents  vorticity  transport  across  the 
boundary  of  the  barrier.  As  a consequence,  barotropic  instability  cannot  be  realized 
in  the  immediate  vicinity  of  the  Tibetan  high  and  the  total  amount  of  energy  extracted 
from  the  larger  scales  by  the  synoptic  scale  within  a fixed  period  of  time  must  be 
greatly  reduced.  Nevertheless,  Krishnamurti  et  al.  (1973a)  were  able  to  simulate  the 
gross  features  of  the  planetary  scale  waves  in  the  upper  troposphere,  such  as  the 
mid-oceanic  troughs  and  the  Mexican  high.  Apparently  once  the  strength,  shape  and 
position  of  the  Tibetan  high  are  fixed  in  the  experiment,  the  interaction  of  the  zonal 
mean  flow  and  this  blocking  high  is  sufficient  to  produce  mechanically  many  other 
features  of  the  circulation. 

In  the  numerical  model  by  Colton  ( 1973),  attention  is  focused  on  the  development 
of  the  synoptic-scale  disturbances  in  the  easterly  Jet  and  in  the  mid-Pacific  trough 
regions  where  the  barotropic  instability  criterion  is  locally  satisfied.  Near  the  easterly 
jet  the  model  produces  disturbances  that  propagate  westward  at  a phase  speed  slightly 
less  than  the  local  speed  of  the  jet  and  continue  to  grow  until  they  are  well  past  the 
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longitude  of  the  jet  maximum.  After  reaching  the  maximum  amplitude  they  begin  io 
decay  fairly  rapidly.  Many  features  of  these  disturbances  resemble  those  analyzed  by 
Krishnaml'rti  and  RtxiERS  ( 1970)  in  the  Indian  Ocean  and  African  region.  In  the 
mid-Pacific  trough  the  model  also  produces  waves  and  vortices  whose  movement  and 
behavior  are  quite  similar  to  those  observed  in  daily  synoptic  charts.  In  both  regions 
the  growing  disturbances  exhibit  a latitudinal  tilt  in  opposition  to  the  lateral  shear  of 
the  quasi-steady  basic  flow,  which  is  characteristic  of  the  barotropically  unstable 
waves. 


Figure  10 

Zonal  average  kineiic  energy  at  20N  as  a function  of  lime  and  wavenumber  as  computed  by  Colton  ( Id73). 
Contour  interval  is  5 m*  sec  Hutched  areas  are  regions  of  energy  maximum. 


In  studying  the  barotropic  energy  exchanges.  Colton  (1973)  also  examined  the 
time  evolution  of  the  computed  kinetic  energy  spectrum.  .An  example  of  this  is  given 
in  Fig.  10  which  shows  the  zonally  averaged  energy  at  20°N  as  a function  of  time  and 
wavenumber.  A noteworthy  feature  is  that  when  the  energy  of  wavenumber  1 is  at  a 
maximum,  the  energy  at  the  short  waves  is  usually  at  a minimum,  and  vice  versa.  This 
suggests  that  energy  exchanges  take  place  continuously  between  the  forced  long  waves 
and  the  transient  disturbances.  The  disturbances  remove  energy  from  the  planetary 
scale  during  the  growing  stage  and  feed  energy  back  into  the  planetary  .scale  during  the 
decaying  stage.  As  a consequence  the  amplitude  of  the  time  mean  planetary-scale 
waves  is  smaller  than  that  produced  by  a linear  model.  Since  no  dissipation  is  provided 
for  the  transient  waves  all  energy  extracted  from  the  forced  motion  must  either 
remain  in  the  transient  waves  or  be  transferred  back  to  the  forced  motion.  Colton 
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(1973)  speculated  that  the  actual  damping  of  the  planetary  scale  by  the  local  baro- 
tropic  instability  may  be  more  effective  if  some  dissipation  mechanism  for  the 
transient  synoptic  waves  is  included. 

Although  Colton's  (1973)  model  lacks  a sufficient  damping  mechanism  and 
therefore  produces  a basic  flow  that  is  somewhat  too  intense  and  has  a phase  shift 
problem,  it  nevertheless  illustrates  the  interesting  phenomenon  of  waves  developing 
and  decaying  in  a spatially  varying  basic  flow  which  contains  locally  unstable  and 
stable  regimes.  To  gain  further  physical  insight  into  this  phenomenon,  especially  with 
respect  to  the  easterly  jet  of  the  summer  monsoon.  Tupaz.  Willia.vis  and  Chang 
(1977)  formulated  a 1-level  numerical  model  using  the  barotropic  vorticity  equation 
to  study  this  problem  from  a linear  instability  point  of  view.  In  their  model  the  basic 
zonal  wind  is  approximated  by  a steady  easterly  Bickley  (hyperbolic-secant-squared) 
Jet  which  increases  slowiy  from  the  inflow  (eastern)  boundary  to  a longitude  in  the 
middle  of  a rectangular  domain,  then  decreases  slowly  downstream  towards  the 
outflow  (western)  boundary.  The  basic  meridional  wind  is  derived  in  such  a way  that 
the  total  basic  flow  is  non-divergent.  Periodic  forcing  is  introduced  on  the  inflow 
boundary  and  a radiation  condition,  whereby  a wave  is  allowed  to  propagate  its 
energy  out  of  the  domain,  is  applied  at  the  outflow  boundary.  The  model  is  integrated 
from  zero  interior  perturbation  streamfunction  until  the  perturbation  fields  vary 
periodically  in  time  at  each  grid  point.  At  this  quasi-steady  state  the  waves  move 
through  the  region  and  grow  or  decay  in  space  only,  in  reaction  to  the  local  stability 
properties  of  the  basic  flow. 

To  compare  the  numerical  results  with  the  ‘local  growth  rate'  concept  of  the 
parallel  flow  theory,  Ti.taz  et  al.  (1977)  developed  a simple  analytical  model  to 
calculate  the  effects  of  the  local  instabilities.  This  is  done  by  considering  an  equation 
that  allows  for  propagation  and  growth  or  decay.  >i/’,  + c{x)\p'^  = n(x)il/'.  where  / is 
time.  .X  the  zonal  coordinate,  and  4/'  the  perturbation  streamfunction.  The  local 
phase  speed  c and  local  growth  rate  n can  be  obtained  from  the  parallel  flow  theory. 
When  a periodic  solution  with  a frequency  to  and  an  amplitude  F{x)  is  introduced,  the 
downstream  variation  of  f with  respect  to  a reference  point  .Vo  may  be  written  as 

Fix)  = F{x„)  exp  [ - ' w t/.v/t(.x)]  exp  [ tt{x)  t/.v/d.x)]. 

Jx„  Jxn 

In  this  equation  the  local  wavelength  is  to/t  and  the  local  spatial  growth  rate  is  o c. 
Figure  1 1 compares  the  steady  state  envelope  for  i/<  of  the  numerical  solution  and  the 
corresponding  Fix)  for  a case  in  which  the  central  speed  of  the  easterly  jet  is  increased 
from  15  m sec" ' at  x = — 10  (XX)  km  to  30  m sec"  ‘ at  x = 0 and  then  reduced  to 
15  m sec"  ‘ at  .X  = 10  (X)0  km.  The  curve  for  Fix)  can  be  shifted  up  or  down  in  the 
diagram  by  changing  x^.  but  its  shape  does  not  change.  Two  Fix)  curves  are  included 
in  Fig.  1 1 , one  for  .Xo  = 280  km  and  the  other  for  Xo  = 4200  km.  both  based  on  the 
local  wavelength  distribution  measured  from  the  numerical  solution.  It  is  evident 
from  this  figure  that  both  the  numerical  solution  and  the  local  growth  rate  of  the 


Vol.  115,  I9''7) 


Some  Theoretical  Problems  of  the  Planetary-Scale  Monsoons 


1103 


figure  1 1 

Companson  of  wave  packet  envelopes  in  a locally  unstable  mean  flow  : ih  . computed  from  a numerical 
model : f jl.Vii  = 280  kmi  and /-^l  v„  = 420<l  kmi.  computed  from  an  analytical  model  based  on  the  parallel 
How  theory  (from  Tli’a/  ct  a/.  1477  j 


parallel  flow  theory  put  the  ma.vimum  amplitude  of  the  perturbations  at  v = -9S00 
km.  well  past  the  jet  maximum.  This  is  in  gotxl  agreement  with  the  results  obtained  by 
Cfii.TON  (1973)  and  ( MANti  and  PfSTiviosn  (I9''7),  and  is  qualitatively  consistent 
with  Krishnavu  Rti  and  RtKii.Rs'  ( I97(|)  observation  that  the  most  active  disturbance 
area  extends  downstream  from  the  jet  maximum  near  the  southern  tip  of  India  to 
central  Africa  (see  Fig.  10.2.  top  panel,  of  Krishsaml  rti.  I97lbi.  This  longitude  of 
maximum  amplitude  corresponds  to  the  neutral  point  of  the  basic  flow,  and  is 
immediately  understandable  from  the  local  growth  rate  concept:  a wave  coming  in 
from  the  inflow  boundary  will  grow  as  it  moves  downstream  in  the  locally  unstable 
region,  until  it  reaches  the  neutral  point  where  it  begins  to  decay.  However,  the 
numerical  solution  has  a larger  maximum  amplitude  than  that  predicted  by  the  local 
growth  rate  of  the  parallel  flow  theory.  The  waves  actually  grow  faster  in  the  unstable 
region  and  decay  faster  in  the  stable  region.  Ti  paz  et  at.  (1977)  evaluated  the  energy 
equation  for  the  numerical  model  and  found  that  the  most  important  term  for 
pyerturbation  energy  production  is  the  Reynold  stress  term  proportional  to  the 
meridional  shear  of  the  basic  zonal  wind.  This  term  is  also  the  only  source  term  for  a 
parallel  flow  model,  and  it  deptends  on  the  phase  structure  of  the  perturbation  field. 
Figure  12  compares  the  local  phase  tilt  of  the  wave  from  the  numerical  model  with 
that  from  the  parallel  flow  theory  at  two  longitudes:  v = 0 km  (jet  maximum)  and 
V = -2SfX)  km  (slightly  downstream  from  jet  maximum).  At  both  places  the  phase 
tilt  from  the  numerical  model  is  significantly  larger  than  that  of  the  parallel  flow  theory. 
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Thus  apparently  the  downstream  variation  of  the  basic  flow  augments  the  phase  tilt 
of  the  waves  which  gives  a larger  growth  rate  in  the  unstable  region.  A similar  situation 
is  found  in  the  stable  region  where  the  dynamic  damping  due  to  the  continuous 
spectrum  effect  is  enhanced  by  a stronger  tilt  in  the  same  direction  of  the  basic  shear  of 
the  easterly  jet.  These  results  indicate  that  the  barotropic  energy  exchanges  associated 
with  a zonally  varying  easterly  jet  which  is  locally  unstable  can  actually  be  more 
prominent  than  those  predicted  by  the  local  growth  rate  from  the  parallel  flow  theory. 
Hence  the  presence  of  the  planetary-scale  asymmetries  during  the  summer  monsoon 
strengthens  the  nonlinear  energy  conversions  resulting  in  significant  activity  of 
synoptic-scale  disturbances. 


Figure  12 

Comparison  of  the  phase  lilts  of  waves  in  a locally  unstable  mean  flow  computed  b>  a numerical  model 
and  the  parallel  flow  theory  longitudes:  tat  x * 0 where  lit  x.  0)  =*  - Ml  m sec  ‘ and 

(b) -x  281)0  km  where  lit  X.  0)  = -27.3mscc‘‘  (From  Tuhaz  <•/ <i/.  1977.) 


Although  Tupaz  ei  at.  ( 1977)  employed  linear  equations  only,  their  results  may 
also  shed  some  light  on  the  eflects  of  the  transient  waves  on  the  basic  flow  . The  waves 
remove  kinetic  energy  from  the  basic  flow  (both  zonal  and  planetary  scale)  and  most 
of  this  energy  is  removed  on  the  downwind  side  of  the  jet  maximum.  Therefore  it  may 
be  expected  that  in  a nonlinear  extension  of  their  model  the  planetary  scale  will 
experience  maximum  damping  west  of  the  easterly  jet  maximum,  which  undoubtedly 
will  alTect  the  amplitude  and  phase  of  the  forced  planetary-scale  waves. 
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4.  Oscillations  of  the  planetary-scale  monsoons 

The  summer  monsoon  circulation  is  known  to  exhibit  many  fluctuations  in  a 
variety  of  time  scales  within  each  season.  The  most  notable  fluctuation  is  the  alter- 
nating active  and  break  periods  of  the  monsoon.  During  the  break  periods  most  of 
the  circulation  characteristics  may  diminish  or  even  reverse  over  the  broad  monsoon 
region.  Less  dramatic  variations  occur  more  frequently  than  the  active-break  cycles, 
although  .ometimes  they  are  limited  to  local  regions.  Figure  13  shows  the  power 
spectrum  of  the  1967  summer  2(X)  mb  streamfunction  over  a region  of  the  Tibetan 
high,  analyzed  by  Krishnaml  rti  et  al.  ( 1973a).  A pronounced  oscillation  periodicity 
is  noted  around  10-13  days  and  a secondary  spectral  peak  occurs  near  3 days.  A 
similar  spectral  distribution  for  the  tropospheric  wind  field  of  the  1962  summer  was 
found  by  M.  Murak.a.mi  (1976),  where  the  main  peaks  occur  at  10-15  days  and  4-5 
days.  Both  of  these  investigators  felt  that  the  shorter  period  oscillation  may  be  a 
manifestation  of  the  synoptic-scale  westward  propagating  monsoon  depressions.  The 
longer  period  oscillation,  on  the  other  hand,  appears  to  be  related  to  the  break  monsoon 
and  other  broad-scale  variations  of  the  monsoon  circulation.  The  10-15  day  periodi- 
city tf  this  oscillation  has  also  been  noted  during  many  other  years.  For  example,  it  is 
reflected  in  thedaily  rainfall  data  of  19b.’.  197]  (.Mi  raka.mi.  1972).  and  1972.  1977 
(Krishna.mi  rti  ana  Bhai.mk.  1976). 

After  detecting  a spectral  peak  near  the  1 3-day  periodicity  in  many  parameters  of 
the  planetary-scale  monsoon  system  of  the  1967  summer.  KRiSHNAviuRriand  Bhaimi; 
( 1976)  calculated  the  phase  relationship  between  the  oscillation  of  these  parameters. 
Their  result  suggests  that  in  a 1 3-day  cycle  the  sequential  strengthening  of  the  mon- 


Figure  I.’ 

Percent  cariance  spectrum  of  itie  mlensily  of  the  Tibetan  High  during  June  August  IR’I . The  dashed  line 
IS  a smoothed  envelope.  (From  Krishnami  rti  «•/ u/..  IV'Ja.i 
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soon  system  including  the  monsoon  trough,  the  cloud  brightness,  the  surface  Mas- 
carene  high,  the  upper  Tibetan  high  and  easterly  jet.  and  the  central  Indian  rainfall 
occur  after  the  dry  and  moist  destabilization  of  the  lower  troposphere.  This  will  then 
be  followed  by  a weakening  of  the  system  in  the  same  sequence  after  dry  and  moist 
stabilization.  On  the  other  hand.  M.  Mlrakami  ( 1976)  used  a composite  technique 
to  study  the  1962  summer  monsoon,  noting  that  it  is  difficult  to  verify  the  statistical 
significance  of  the  cross-spectral  quantities  of  this  period  range  derived  from  a 
3-month  time  series  of  daily  data.  He  found  that  the  intensification  and  weakening  of 
convection  and  other  low-  and  upper-level  winds  are  closely  in  phase  for  the  15-day 
cycic  during  this  season.  In  any  case,  there  is  strong  evidence  that  many  components  of 
the  planetary-scale  monsoon  are  involved  in  this  oscillation. 

The  observed  10-15  day  oscillation  raises  an  interesting  question  concerning  its 
physical  mechanism.  Krishnamlrti  and  Bhal.vii;  (1976)  proposed  that  it  may  be  a 
result  of  a cloud-radiative  feedback  mechanism  in  which  the  enhanced  convection 
causes  increased  cloud  cover  which  in  turn  reduces  the  incoming  short-wave  radiation 
and  stabilizes  the  lower  troposphere.  This  is  then  followed  by  decreased  convection 
and  cloudiness,  and  therefore  increased  short-wave  radiative  heating  at  the  surface 
which  renews  the  active  convection.  Another  possible  oscillation  mechanism  is 
revealed  by  a simple,  domain-averaged,  atmosphere-occan  interacting  model  by 
Webster  and  Lau  (1977).  In  their  model  the  large  heat  capacity  of  the  ocean  causes 
it  to  respond  slowly  to  the  solar  heating  while  the  atmosphere  temperature  reaches  a 
ma.\imum  due  to  heating  at  the  land  surface.  After  a period  of  time  the  ocean  tem- 
perature peaks  and  this  feeds  back  to  the  atmosphere  through  convection  resulting  in 
another  temperature  maximum.  Although  this  process  alone  is  unlikely  to  produce 
subsequent  maxima  without  the  participation  of  other  mechanisms,  it  may  still  be 
relevant  to  some  of  the  active-break  cycles  because  Ra.ma.mijRTHy's  (19h9)  survey 
' (quoted  in  Krishnamurti  and  Bhalvie.  1976)  indicates  that  there  are  only  one  or  two 
major  breaks  during  most  of  the  summer  monsoon  seasons.  Other  possible  mechan- 
isms for  the  planetary-scale  oscillation  may  include  interactions  between  hemispheres 
and  between  stationary  and  propagating  long  waves  as  mentioned  by  Krishnamurti 
and  Bhalme  (1976).  and  nonlinear  vacillations  between  different  scales.  The  latter 
mechanism  produces  a 5-day  oscillation  of  the  planetary  scale  in  Colton's  (1973) 
barotropic  model,  which  is  purely  dynamical  and  does  not  provide  dissipation 
proces.ses  for  the  synoptic  scale.  If  thermodynamics  and  synoptic  damping  are 
included,  the  interaction  would  become  more  complex  and  the  oscillation  period 
could  be  modified. 

The  planetary-scale  circulation  during  the  winter  monsoon  also  exhibits  significant 
variations,  largely  due  to  the  influence  of  the  cold  air  outbreaks  from  the  Asian 
continent.  These  outbreaks  often  interfere  with  the  main  winter  convective  region  in 
the  maritime  continent  and  cause  widespread  weather  changes  in  the  South  China 
Sea  and  .Vlalaysia-lndonesia  region,  including  the  development  of  strong  synoptic- 
scale  disturbance  activity.  Since  this  convective  region  provides  the  major  energy 
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source  for  the  thermally  direct,  planetary-scale  east-west  overturning  and  Hadley 
cells  (Krishna.mlrti  et  al.  1973b).  its  variations  directly  influence  the  circulations  in 
the  entire  tropics  and  southern  mid-latitudes  during  winter.  These  planetary-scale 
zonal  and  lateral  teleconnections  are  reflected  in  the  analyses  by  Muraka.mi  and 
Unninayar  (1977)  and  Krishnamurti  ft  al.  ( 1973b)  despite  their  limited  data  base. 
The  initiation  and  progress  of  theoretical  studies  in  this  area  would  depend  heavily  on 
the  expansion  of  our  present  observational  knowledge  of  this  problem. 


5.  Concluding  remarks 

The  three  problems  reviewed  here  are  inter-related  to  each  other.  For  example, 
the  localized  barotropic  instability  is  a reflection  of  the  strong  nonlinear  energy 
transfer  which  appears  to  be  quite  relevant  to  the  damping  of  the  planetary-scale 
flow  during  the  summer  monsoon.  Also  Chang  and  Pentimonti's  ( 1977)  experiment 
indicates  that  these  energy  conversions  are  more  pronounced  when  the  planetary 
scale  contains  an  appreciable  oscillation.  The  enhanced  upper  tropospheric  in- 
stability in  summer  due  to  both  time  and  spatial  variations  of  the  planetary  scale  also 
has  important  ramifications  for  the  synoptic-scale  disturbances.  At  least  in  the  mid- 
oceanic  trough  regions,  the  upper-level  disturbances  may  develop  in  the  vertical  and 
could  be  related  to  some  of  the  weather-producing  easterly  waves  and  tropica) 
cyclones.  A ditferent  type  of  problem  arises  for  the  winter  monsoon  during  which  the 
upper  tropospheric  synoptic-scale  waves  and  zonal  flow  feed  energy  to  the  planetary 
scale  (Krishna.mlrti  ct  al.  I973bi.  This  on  the  one  hand  removes  the  nonlinear 
proces.ses  as  a possible  damping  mechanism,  and  on  the  other  hand  requires  an 
explanation  for  the  energy  source  of  the  synoptic  scale.  Clearly  much  more  elfort  is 
needed  to  answer  these  questions,  whose  importance  lies  in  the  fact  that  a thorough 
understanding  of  the  planetary  scale  is  crucial  for  the  proper  simulation  of  the  entire 
multiple-scale  monsoon  circulation. 
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